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PREFACE 


During the International Association of Geomagnetism and Aeronomy (IAGA) General Scientific Assembly in Madrid, 1-12 
September 1969, Commission V (Solar-Magnetosphere Relations) held a one-and-a-half day Scientific Session on Models of the 
Earth’s Radiation Environment. This was organized by the undersigned at the suggestion of the Chairman of Commission V, Dr. 

J. G. Roederer. In addition to quantitative models of the radiation environment, long-term time variations and average behavior of 
magnetospheric particles were emphasized. The following papers were presented by the authors whose names are In italics. 

Imhof W. L . : “The Low Altitude Radiation Environment” (invited), 

Bostrom , C. 0 ., Beall, D. S., and Armstrong, J. C.: “Time Variations in the Inner Radiation Zone,” 

Rosen , A. and Sanders, N. L.: “The Temporal Behavior of Inner Zone Electrons During 1965, 1966 and 1967,” 

Heckman , H. H. and Nakano, G. H.: “Long-Term Behavior of Energetic Inner Belt Protons,” 

Filz, R. C. and Bellew, W. F. \ “Comparison of Trapped Radiation Models with Nuclear Emulsion Measurements at Low 
Altitudes,” 

Rothwell, R: “Energetic Particles and Plasma Inside the Earth’s Plasmasphere” (invited), 

Frank , L. : “Spatial Distributions of Low Energy Protons and Electrons in the Earth’s Magnetosphere” (invited), 

Davis, L. R.: “Spatial and Temporal Characteristics of 100 KeV to 1700 KeV Outerbelt Protons” (invited), 

Webb, V. H, Rothwell, P. L., and Katz, L.: “Quiescent and Disturbed Proton and Electron Distributions,” 

Sanders, N. L. and Rosen , A. : “Study of the Response of E > 0.3 MeV and E > 0.5 MeV Electrons in the Inner Zone to 
Magnetic Activity,” 

Vette, J. I. and Singley, G. W.: “An Outer Zone Electron Model Environment” (invited), 

Paulikas, G. A. and Blake, J. B. : “The Particle Environment at Synchronous Altitude” (invited), 

White, R. S Wilson, J., Martina, F., and Stevens, J. R.: “Proton Energy Distributions from 0.065 to 3.6 MeV at 6.6 Earth 
Radii,” 

Pruett, R. G, and Lucero, A. B.: “The Impact of the Charged Particle Space Radiation Environment on Spacecraft Design and 
Mission Analysis,” 

Nishida, A., Lyon, E. F., and Ness, N. F. : “Plasma Sheet at Lunar Distance,” 

Anderson, K. A.: “High Energy Particles in the Earth’s Magnetotail ” (invited), 
and 

Masley, A. J.: “Solar Cosmic Ray Activity Near Sunspot Maximum.” 

A number of these papers were solicited for publication in this series of Models of the Trapped Radiation Environment be- 
cause they contained an excellent description of the long-term behavior of several components of the radiation belts. The four 
papers presented here are the result of that effort. I wish to thank all the contributors to the sessions in Madrid and especially the 
authors who have provided their manuscripts for publication herein. The Outer Zone Electron Model Environment by G. W. 
Singley and myself will form the main body of the next volume in this series. 

James I. Vette 

National Space Science Data Center 
Goddard Space Flight Center 
Greenbelt, Maryland 
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THE LOW-ALTITUDE RADIATION ENVIRONMENT 


W. L. Imhof 

Physical Sciences Laboratory 
Lockheed Palo Alto Research Laboratory 
3251 Hanover Street 
Palo Alto, California 94304 




ABSTRACT 


Observations of electrons and protons at low altitudes have provided valuable information relat- 
ing to the source and loss mechanisms of particles trapped in both the inner and outer radiation belts. 
In the various latitude regions, the radiation environment is quite different in character. At the lower 
edge of the inner radiation belt the flux gradients of both electrons and protons are typically quite 
steep and the lifetimes rather short. Accordingly, this region of space has been important for evaluating 
the cosmic-ray albedo-neutron decay source of trapped protons. Observations of redistributions of 
trapped electrons associated with magnetic storms and measurements of the long-term decay of the 
Starfish radiation belt have indicated the occurrence of L-shell diffusion. Throughout the inner radia- 
tion belt the fluxes of trapped electrons having minimum altitudes below sea level are seen to vary with 
longitude, indicating the continual action of precipitation processes. In the slot region and the outer 
radiation belt the rates of buildup with longitude drift of the quasi-trapped electrons have been found 
to be consistent with the assumption that the dominant loss mechanism for trapped electrons is pitch- 
angle diffusion. More detailed information on the precipitation mechanisms has been obtained from 
direct observation of electrons having pitch angles within the loss cone. In the polar-cap regions, low- 
altitude observations of the incident particle fluxes have revealed a great complexity in the parameters 
associated with energetic electrons, protons, and alpha particles. The foregoing topics are reviewed 
with emphasis on the long-term average behavior. 


INTRODUCTION 


The low-altitude radiation environment encompasses many regions: the inner and outer radiation belts, the auroral zones, and 
the polar caps. In any of these regions, at heights up to a few hundred kilometers, the atmosphere strongly influences trapped and 
precipitating particles. As a result, the spatial flux gradients are often very steep and the trapping lifetimes quite short. These fea- 
tures, when measured in detail, can provide unique clues to the important source and loss mechanisms for particles trapped in the 
earth’s radiation belts. The purpose of this paper is to present some of the more important low-altitude measurements and to at- 
tempt to summarize the environment and related geophysical characteristics obtained from these observations. In order to cover the 
subject with any reasonable degree of completeness, the discussion is restricted primarily to those characteristics related to the higher 
energy (^40 keV) electrons and (^40 MeV) protons, for which more complete data are available than for the lower energy fluxes. 
Even under such restrictions a complete review of all experiments at low altitudes is clearly not possible within the scope of this paper. 


INNER RADIATION BELT 

The distinction between the earth’s inner and outer radiation belts (with a slot region in between) is not precisely defined. 

The change in character with position depends strongly upon the energy and type of particle. A definite slot region exists for MeV 
electrons but it tends to disappear at lower energies. High-energy protons, on the other hand, exist only in the inner radiation belt. 
These general features for both electrons and protons are summarized in Figure 1, where the equatorial fluxes are plotted as a func- 
tion of Mcll wain’s (1961) ^-parameter. In order to discuss the low-altitude portion of each of these regions in some logical order, 
we start with the inner belt, progress outward through the outer belt, and finally consider the polar-cap regions. In starting with 
the inner belt, much of the initial attention is devoted to what is perhaps the most distinctive low-altitude region of the radiation 
belts, the South Atlantic anomaly. Here the earth’s magnetic field is distorted so that trapped particles dip to the lowest altitude 
in their longitude drifts around the world. 

Trapped Protons 

Much of our knowledge of inner-zone protons has come from measurements performed at relatively low altitudes. Some of 
the earliest and most significant of these measurements were made with nuclear emulsions carried by recoverable rockets or 
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satellites. For example, the first identification of protons in the inner belt resulted from the emulsion experiment of Freden and 
White (1959). The energy spectra obtained from this and similar flights are shown in Figure 2a. The shape, which agrees rather 
well with subsequent measurements is consistent with the expectation for injection of trapped protons by the GRAND (cosmic-ray 
albedo-neutron decay) process. The agreement in spectral shape becomes less valid at lower proton energies, particularly for 

T-shells above about 1.4. 

At low altitudes, pronounced changes with solar-cycle period in the energy spectrum of high-energy protons have been pre- 
dicted by Blanchard and Hess (1964). Thus, one should expect the low-altitude trapped proton environment to be very complex 
and interesting in regard to spectral and intensity variations. However, as shown in Figure 2b (Heckman and Nakano, 1969), ob- 
servations have indicated a lack of pronounced time variations (through June 1966) in the proton spectrum, perhaps due to the 
small atmospheric changes which have taken place during this period. 

One of the more dramatic characteristics of the trapped proton environment at low altitudes is the existence of an east-west 
asymmetry, shown schematically in Figure 3. The asymmetry, first predicted by Lenchek and Singer (1962), arises because (at the 
point of observation) protons arriving from the east have guiding centers at lower altitudes than those arriving from the west. Con- 
clusive observations of this asymmetry were first reported by Heckman and Nakano (1963). In their experiment protons above 57 
MeV were recorded in nuclear emulsions placed on a low-altitude oriented satellite. The angular distributions of protons stopping 
in an emulsion stack displayed unequal numbers of protons entering from the east and west. This effect has also been observed 
with emulsions by Filz and Holeman (1965); and with counters by Garmire (1963), Galperin and Temny (1965), and Reagan and 
Imhof (1969a). Overall, these various observations show general agreement with the assumption that the dominant loss mecha- 
nism is atmospheric interactions. However, in detail certain deviations are observed. Heckman and Nakano (1969) have found 
that the asymmetries differ significantly from calculations for the Harris and Priester model atmosphere. For such comparisons 
with theory the asymmetry measurements provide information similar to that derived from omnidirectional flux contours, the 
latter being used more widely. The east-west asymmetry has yet to be measured in great detail as a function of position in space. 
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cosmic-ray albedo-neutron decay. (Data from F reden 
and White, 1960.) 


Figure 2— Energy spectra of inner-zone protons. 

The absolute intensities of protons trapped in the inner radiation belt have been compared by several authors (Haerendel, 
1964; Hess and Killeen, 1966; Dragt et al., 1966; and Imhof and Reagan, 1969a) with calculations of the GRAND source strength 
and atmospheric losses. The calculated GRAND source strength appears to be too weak throughout the inner belt by factors rang- 
ing from 6.5 to 50. It has been argued that the ideal positions for comparing the measured fluxes with the GRAND expectations 
are at lower altitudes where any diffusion effects should be small. Accordingly, Dragt et al. ( 1 966) compared their calculated pro- 
files of flux versus B at L ~ 1.4 with low-altitude emulsion measurements performed by Filz and Holeman (1965). As is shown In 
Figure 4, a discrepancy of a factor of 50 between CRAND and the measured fluxes resulted. In an attempt to reduce further any 
pitch-angle and/or T -shell diffusion contributions, Imhof and Reagan (1969a) have considered the lower edge of the Inner radia- 
tion belt. At L < 1 . 1 5 the slopes of the spatial flux profiles were found to be in reasonable agreement with the predictions of an 
injection rate nearly independent of position over this narrow region of space and a loss rate due to atmospheric interactions, as 
shown in Figure 5. However, the measured fluxes were higher by about a factor of 20 than those expected for CRAND. 

Effects of atmospheric variations on trapped-proton fluxes have been observed by Nakano and Heckman (1968). Between 
November 1962 and June 1966, the omnidirectional proton flux at 63 MeV remained constant to within ±7.6%, as shown in Fig- 
ure 6. Then, between August 1966 and November 1967, the proton flux decreased by about a factor of two. A drop in flux at 
this time in the solar cycle was predicted by Blanchard and Hess (1964), but a much larger decrease was expected. However, this 
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Figure 3— Schematic representation of the east-west asymmetry. 
Due to the atmospheric density gradient, the fluxes of protons 
traveling toward the west are lower than those traveling toward 

the east (Heckman and Jakano, 1965). 

apparent discrepancy may not be significant since the actual 
solar-cycle variations of the 10.7-cm flux have been consider- 
ably less than those used in the calculations of Blanchard and 
Hess. Equatorial fluxes of protons were measured over the 
range 1.15 <L< 1.21 by Achtermann et al. (1969). Their 
measured profiles of flux versus L displayed variations at dif- 
ferent proton energies, and hence different lifetimes, which 
they attribute to changes in the atmosphere shortly before the 
flight on June 17, 1967. 



Figure 4— Comparison of the proton fluxes measured by Filz 
and Hofeman with the calculated fluxes for cosmic-ray albedo- 
neutron decay (Dragt et al., 1966). 


Many of the fundamental questions relating to the sources 
of protons in the inner belt and possible redistribution mecha- 
nisms are more properly considered elsewhere. Here some of the effects unique to the low-altitude regions have been outlined. 
Several classes of experiments have shown the dominance of atmospheric effects. The east-west asymmetry, the spatial flux pro- 
files, and the variations of the trapped-proton fluxes with periods in the solar cycle are generally consistent with the assumption 
that the loss rate is determined by atmospheric collisions. However, significant deviations have been observed. Perhaps more seri- 
ous discrepancies exist in the absolute fluxes; even at the lower edge of the inner radiation belt the intensities of high-energy pro- 
tons are much larger than predicted for the GRAND source mechanism. 


Trapped Electrons 

Very few measurements of electrons in the inner radiation belt were made before the artificial injection of electrons by the 
Starfish nuclear detonation on July 9, 1962. In particular, the energy spectrum was not well established. The available pre-Starfish 
electron spectrum measurements, as illustrated in Figure 7, display a large degree of inconsistency. Intensity profiles are available 
only for limited spatial regions. It is clear that our knowledge of the naturally trapped electrons in the inner belt will remain mea- 
ger until the Starfish electrons no longer dominate. Fortunately, the emergence of the natural electron belt may be just around the 
comer; recent experiments have indicated the appearance of naturally trapped electrons in the inner belt (Marbach et al., 1968); 
Pfitzer and Winckier, 1968; and Vampola, 1968). But at present our knowledge of the behavior of electrons in the inner belt is 
derived mostly from measurements of Starfish electrons. 
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Figure 5~Equatoriai omnidirectional fluxes plotted as a func- 
tion of L. The points represent measured values with 
150 km normalized to the equator on the basis of the calcu- 
lated flux variations with B. The curve represents the fluxes 
calculated for the GRAND source strength arbitrarily increased 
by a factor of 20 (Imhof and Reagan, 1969a). 


The sudden injection of many electrons from the Starfish 
detonation presented a unique opportunity to measure trap- 
ping lifetimes. The early lifetime observations performed on 
Injun 1 (Van Allen, 1963) and Explorer 15 (Mellwain, 1963) 
have been compared with those calculated by Walt (1964) for 
losses of electrons due to collisions with atmospheric constitu- 
ents. For L ^ 1.25, the measured lifetimes were considerably 
shorter than the theoretical values, suggesting the dominance 
of other loss processes. However, as is shown In Figure 8, for 
L ^ 1.25 there was good agreement between theory and exper- 
iment, indicating that atmospheric scattering is the dominant 
loss mechanism in this region of space. The lifetime measure- 
ments were extended down to an A-value of 1. 145 by Imhof 
and Smith (1966). They observed the decay of electron inten- 
sities following two different energy-selective redistributions of 
electrons at 1.34 and 0.74 MeV. Each of the two groups of 
electrons decayed at a rate that was in agreement with the at- 
mospheric scattering calculations of Walt. 

During a somewhat later period (October 1 963 through 
December 1965) the decay of the Starfish radiation belt at low 
altitudes was followed in detail with instruments aboard the 
satellite 1963 38C (Beall et a!., 1967). The long-term averages 
displayed a systematic dependence on energy, time, and posi- 
tion, which for L » 1.25 was in general agreement with Walt’s 
calculations at relatively early times when a 5-depen dence in 
decay rate was expected. However, the observed decrease in 
lifetime with increasing B was relatively weak. Despite the 
overall average agreement with atmospheric scattering calcula- 
tions, the observed decay was somewhat irregular; significant 
short-term variations in the decay constant were observed. 

Electrons trapped on low L-shells have been observed to 
undergo significant variations in intensity. Measurements from 
the Pegasus 1 satellite (Rosen et al., 1968) in 1965 revealed an 
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Figure 6— Time behavior of the proton flux between September 1962 and November 
1 967. The fractional deviations are plotted as a function of time for 220 < h . < 375 

min 

km and 375 </? min < 455 km (Nakano and Heckman, 1968). 
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Figure 7— Pre-Starfish electron energy spectra (Imhof et a!., 1963). 

enhancement of 0.5-MeV electron fluxes at Z-shells less than 1.2. This enhancement was related to, but not clearly correlated 
with, magnetic activity. Measurements on a polar-orbiting satellite with an instrument having high-energy resolution and sensitivity 
have shown, on two different occasions, the appearance of nearly monoenergetic groups of electrons (Imhof and Smith, 1965b and 
1966) on low £ -shells. One of the sets of observations is illustrated in Figure 9. Data from the satellite measurements and from 
ground magnetogram records suggest that the observed onset of nearly monoenergetic electrons was the result of a redistribution 
of electrons trapped on somewhat higher £ -shells (Imhof et al., 1966). The responsible mechanism has been explored by Cladis 
(1966). He has suggested that the electrons were accelerated by fluctuations in the intensity of the equatorial electrojet that had 
periods comparable to the azimuthal drift periods of the electrons. The energies imparted to the electrons were small but suffi- 
cient to account for the observed redistributions. For this acceleration mechanism the energy spectra computed by Cladis were in 
good agreement with the measurements. 

It is evident that in the absence of any sources of new electrons the lower £-shells should eventually become completely 
vacated due to atmospheric scattering. However, measurements (Imhof et al., 1967 and Rosen et al., 1968) have indicated that 
three years after Starfish, electrons were still trapped on £-shells ^ 1.18. The persistence of trapped electrons at the lower edge of 
the inner radiation belt has been interpreted by Newkirk and Walt (1968) in terms of diffusion across £-shells. Empirical values 
for the diffusion coefficients at low £-shells have been obtained by Newkirk and Walt (1968) and by Farley (1969a and b). Both 
of these calculations were based on the measured radial-flux profiles and long-term decay rates. The values of the diffusion coeffi- 
cients were found to increase monotonically with decreasing £, a trend which is opposite to the variation of diffusion rates higher 
in the magnetosphere. The physical interpretation of this unique characteristic of the very low £-shell region is still an open 
question. 

From the collection of observations at low £-shells, it appears that the trapping lifetimes at £ ^ 1.25 are consistent with the 
assumption that the dominant loss mechanism is atmospheric scattering. On much higher £-shells, clearly other loss mechanisms 
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Figure 8— Electron decay rate as a function of L. The time constant is the time required 
for the electron flux to decrease by a factor of Me. The theoretical curve for L = 1 .185 
was calculated for an initial fission energy spectrum; the theoretical curve near L = 1.15, 
for a nearly rnonoenergetic group of electrons at 1 .3 MeV (Walt, 1966). 

are dominant. On very low L-shells the behavior of the electrons is rather complex. Atmospheric scattering appears to dominate 
at times but other perturbing effects play an important role. Significant redistributions have been reported at times of magnetic 
disturbances. Overall, there appears to be a significant radial movement of electrons from one i-shell to another. 

One of the outstanding features unique to the low-altitude environment is the existence of electrons with mirror points at 
positions where the longitude trace dips to an altitude below sea level. Such electrons may appear to be trapped, but are doomed 
to precipitate into the atmosphere as a result of their longitude drift. Accordingly, these fluxes of quasi-trapped electrons mirror- 
ing at a given B-, Z -value are seen to display pronounced variations with longitude. The combination of the flux buildup with 
longitude drift and the subsequent dumping is often referred to as the wind shield- wiper effect. This phenomenon has been re- 
ported over a very wide range of Z-values (Paulikas and Freden, 1964; Imhof and Smith, 1965a; Vemov et aL, 1965; Williams and 
Kohl, 1965; Blake et aL, 1966; Rosen et al, 1968; and Imhof, 1968). As an illustrative example, the fluxes of electrons above 50 
keV measured by Vemov et al. (1965) on Cosmos 17 in mid- 1963 are shown in Figure 10. On various Z-shells the intensities are 
plotted as a function of h min - As might be expected for electrons drifting eastward, the fluxes are greater to the west of the Brazil- 
ian anomaly than to the east. The appearance of particles at low values of h min appears to be more prevalent for electrons than pro- 
tons; however, Vemov et al. (1966) have reported a longitude dependence of protons on magnetic shells L = 1.2 to 1.4 in the 
region of the South Atlantic anomaly. 

It can be seen from Figure 10 that in the inner belt the fluxes fall off very rapidly with decreasing h min in the vicinity of 
h m in ~ 100 km. Whereas the electrons at h mhl ~ 100 km may represent an appreciable fraction of those being lost by atmospheric 
scattering from the belt, the intensities at h min < 0 km are considerably lower. It has been shown (Imhof and Smith, 1965a) that 
in late 1963 on a given Z-shell in the inner belt, the buildup and subsequent precipitation of electrons at h min < 0 km can be sus- 
tained by an influx considerably less than the total loss rate of trapped electrons from that shell. Consequently, the fluxes appear- 
ing at high ^-values may be associated with a small effect which is perhaps not representative of the primary loss mechanism(s). 
Atmospheric scattering calculations by Roederer et al. (1967) for low Z-shells have shown that multiple coulomb scattering can 
account for the electron fluxes observed at ^ min ^ 0 km, but that an additional mechanism is needed to explain the fluxes observed 
at h. < 0 km. 
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Figure 9—Output of an electron spectrometer on different 1-shells for one pass through the South Atlantic anomaly 

(Imhof and Smith, 1965b). 


Various processes have been suggested to account for the electrons observed at minimum altitudes below sea level, but for any 
of these mechanisms detailed comparisons between theory and experiment are still lacking. On the basis of the intensities meas- 
ured in September 1962 and in October to November 1963, it has been suggested (Imhof and Smith, 1965a) that inner-belt elec- 
trons at & min < 0 km originate from the trapped-radiation belt. An important redistribution process was postulated to be single 
large-angle coulomb scatterings in the atmosphere. Such scatterings surely occur, but more extensive calculations and comparisons 
with the measured intensities and longitude variations are needed to firmly establish the importance of this process. Other non- 
atmospheric precipitation mechanisms have been suggested, particularly to account for the high-i? fluxes in the slot region and 
outer belt where the atmosphere clearly cannot play an important role. The scattering of trapped electrons by whistlers was first 
investigated by Dungey (1963) and by Cornwall (1964). Cornwall calculated that the conditions for a resonance interaction lead 
to a pronounced dependence on I, the maximum interactions occurring near L = 1.2 and L = 2. Subsequently, Roberts (1968) 
has performed calcuiations to support the suggestion of a dual mechanism to account for the precipitation of trapped electrons. 

He proposed that bounce-resonant scattering accounts for the pitch-angle diffusion of electrons mirroring near the equator, 
whereas those mirroring at lower altitudes are perturbed by cyclotron-resonant scattering. Insufficient data, particularly on the 
magnetic and electric field fluctuations in the magnetosphere are presently available to evaluate definitively this dual mechanism. 


OUTER RADIATION BELT 


Many satellite measurements of electrons in the outer radiation belt have been made at low altitudes. Perhaps one of the 
most notable sets of measurements in regard to long-term coverage has been performed from 1963 38C (October 1963 to present). 
In this paper, processes such as precipitation phenomena that are more characteristic of the low-altitude regions are emphasized. 
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ALTITUDE, h m j n (km) 

Figure 10— Intensities on various magnetic shells plotted as a 
function of minimum altitude of the mirror points. The open 
circles represent measurements performed at positions east of 
the South Atlantic anomaly (Vernov et al., 1965). 

electrons are injected into the depleted regions at all longitudes, 
mission through the South Atlantic anomaly as a function of h n 
h . > 350 km, about 10% for h ■ = 100 + 250 km, and 0% foi 


As already mentioned in the discussion on the inner radi- 
ation belt, quasi- trapped electrons doomed to precipitate upon 
drift in longitude have been regularly observed throughout the 
radiation belts. Paulikas et al. ( 1 966) have noted that such 
electrons appear in three broad peaks in L ~~ 1.5, ~ 2.5, and 4 
to 7. In the slot region and the outer belt, the character of the 
electron precipitation is in distinct contrast to that measured 
in the inner belt. On low L-shells significant fluxes of elec- 
trons at low altitudes are generally limited to the immediate 
vicinity of the South Atlantic anomaly. On the other hand, in 
the outer belt appreciable fluxes of precipitating electrons are 
observed to exist on a worldwide basis (Paulikas and Freden, 
1964 and Imhof, 1968), as is shown in Figure 1 1. in the inner 
belt the energy spectrum of electrons with h m - n < 0 km lias 
been observed to be much softer than that of the more stably 
trapped electrons. However, in the outer belt the spectra dis- 
play only minor variations, with h min being somewhat softer 
for the precipitating electrons (Imhof, 1 968). In the inner belt 
the intensities of the quasi-trapped electrons at any given re- 
gion of space have been observed to be reasonably steady. 
Williams and Kohl (1965) have noted in the inner-belt and slot 
regions that a relatively stable longitudinal intensity distribu- 
tion exists during periods of magnetic quiet. However, signifi- 
cant fluctuations have been observed in the slot and outer-belt 
regions (Paulikas and Freden, 1964 and Imhof et al., 1969). 
Despite the scatter in fluxes of precipitating electrons, both 
the findings of Williams and Kohl (1965) and of Imhof (1968) 
have indicated a general flux buildup with eastward drift in 
longitude. 

To date, studies of the flux variations with longitude have 
been quite limited in number and scope. From a series of flux 
measurements at several longitudes in the northern hemisphere, 
Williams and Kohl (1965) concluded that for 1.45 < L < 2.8 
They also obtained measurements of the fractional electron trans- 
n . The transmissions were found to be 1 00% for electrons with 
h. < 100 km. 


The buildup of electron flux with longitude drift has been investigated with data taken in 1963 on a polar-orbiting satellite 
that had complete worldwide coverage through use of an onboard tape recorder (Imhof, 1 968). The longitude variations of the 
intensities and spectra were found to be consistent with a continual buildup and, in the outer-belt regions, with the rapid achieve- 
ment of a state of near equilibrium. The rates of flux buildup were compared with calculations of pitch-angle diffusion, yielding 
empirical values for the diffusion coefficients. These coefficients were found to vary by only a factor of four over a rather wide 
range of magnetic latitudes covering the outer radiation belt and slot regions. The diffusion coefficients were also found to be con- 
sistent with the assumption that the dominant loss mechanism for trapped electrons is pitch-angle diffusion. Pitch-angle diffusion 
coefficients for electrons in the outer belt have also been obtained (Theodoridis and Paolini, 1967) by measuring the pitch-angle 
distribution near and within the loss cone. The diffusion coefficients derived by this method were considerably larger than those 
obtained by Imhof from the buildup of the trapped flux with longitude drift. 


The observations of electrons in the loss cone are most conveniently performed at low altitudes since here the loss cone is 
relatively large. Some of the earliest detailed studies of the precipitation of electrons at high latitudes were performed with instru- 
ments onboard the satellites Injun 1 and Injun 3 by O’Brien (1964). He studied in detail the variations in the intensities of the 
trapped and precipitating electrons above 40 keV. During precipitation events, the angular distributions tended to approach isot- 
ropy over the upper hemisphere, as is illustrated in Figure 12. Strong energy dependences were observed. O’Brien found that elec- 
trons > 1 MeV can be trapped and unperturbed on a field line along which fresh electrons are being precipitated with energies 
above 40 keV. 
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Figure 11 — Locus of peaks of electron flux above 0.9 MeV for September 1962 as observed at altitudes of 300 to 
690 km. Measurements near the center of the South Atlantic anomaly are not included {Paulikas and Freden, 1964). 

Further evidence for the tendency to favor the precipitation of low-energy electrons has come from the observations ofMozer 
(1 968). From high-energy (> 65 keV) electron measurements on rockets flown into auroras at times near local midnight, Mozer 
found that the variations were generally greatest for electrons with energies between 65 and 200 keV. However, fluctuations occa- 
sionally occurred exclusively in the electrons above 400 keV. During five well-defined magnetic storms, Williams et al. (1968) 
compared simultaneous flux measurements of electrons at low altitudes and in near-equatorial regions, as is shown in Figure 13. 

The behavior of electrons of about 300-keV, mirroring near the equator and at low altitudes is very similar; the >300-keV electron 
intensities on a given L-shell generally come quickly to equilibrium. On the other hand, for electrons > 1-MeV such a rapid attain- 
ment of equilibrium is not achieved, suggesting that the mechanisms of major importance in pitch-angle diffusion apparently act 
more effectively for electrons >300-keV. 

Energy dependences of the overall loss rates of trapped electrons can be obtained from the measured time variations in the 
intensities of the trapped population. Vernov et al. (1969) have concluded from all available measurements that in the 40- to 500- 
keV energy range the lifetime depends little on energy, but that above 500 keV it increases rapidly with increasing energy. From 
all of these classes of experiments one concludes that on the average low-energy electrons precipitate at a somewhat greater rate. 

The foregoing evidence for some energy selectivity in the precipitation of electrons leads to the interesting considerations of 
the existence at high latitudes of nearly monoenergetic groups of electrons, as observed at the lower edge of the inner radiation 
belt. Observations show that the lower energy auroral electrons do display monoenergetic groups. Mcllwain (1960) first suggested 
the presence of a monoenergetic peak at about 6 keV in interpreting some of his auroral zone measurements. More recently, nearly 
monoenergetic groups of electrons at a few keV have been observed by Evans (1968), Albert (1967a and b), and by Westerlund 
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Figure 12— Fluxes of trapped and precipitating electrons plotted as a function of 
invariant latitude for one pass. The outputs of a photometer onboard the satellite are 
also shown (O'Brien, 1964). 


(1969). High-altitude observations in the slot region of spectra peaked at energies of several hundred keV have been reported by 
West (1969). 

Overall, the measurements of high-energy electrons at high latitudes have been rather limited in spatial and temporal coverage, 
sensitivity, and in energy and angular resolution. These deficiencies can be attributed to the need for fast-time resolution and to 
the relatively low flux levels present. Greater intensities are generally present at lower energies, but until recently the much greater 
experimental difficulties associated with measuring low-energy particles have severely limited the coverage. Now with significant 
advances in experimental techniques, this situation has changed somewhat, but there still exist very few systematic studies of the 
overall morphology as measured from satellites. These measurements are summarized in recent review articles (Sharp and Johnson, 
1968 and Johnson and Sharp, 1969). 

Some of the most interesting and important features of the low-altitude radiation environment involve the outer trapping 
boundary. Several studies of the high-latitude boundary have revealed the complexities of this region. The boundary is not only 
dependent upon various parameters such as the local magnetic time and the magnetic activity, but at any one time its location is 
strongly dependent upon the precise definition being used. Perhaps the simplest definition of a boundary in terms of data analysis 
is the position at which a given detector reaches a specified counting rate. Such a prescription is convenient for long-term studies 
with a given set of instruments but is clearly dependent upon the response of the detector and does not necessarily correspond to a 
sudden change in the physical character of the radiation environment. 

Williams and Mead (1965) have studied diurnal variations of the electron-trapping boundary using absolute electron intensities 
as a latitude criterion. The diurnal variations of the boundaries for electrons above 280 keV and above 1.2 MeV were found to be 
similar and consistent with the drift of the particles through a distorted magnetosphere. By comparing the calculations with 
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experiment they obtained an empirical nightside magnetic field configuration. From comparisons with similar measurements of 
electrons >40-keV by Frank et a!. (1964), they found that the midnight boundaries were nearly independent of electron energy 
but that the dayside boundary for the low-energy electrons extended to significantly higher latitudes than was the case for the 
higher energy electrons. Williams and Mead interpreted the larger diurnal shifts for electrons >40-keV to imply a daytime source 
mechanism. More recently, Taylor (1966) has accounted for the differences in daytime latitudes of the boundaries >4G-keV and 
>280-keV by invoking magnetospheric electric fields. 

With data acquired on the Alouette 2 satellite the Canadian group has examined a variety of high-latitude boundaries (Burrows 
and McDiarmid, 1968 and McDiarmid and Burrows, 1968). Each of these boundaries is defined in terms of the absolute counting 
rates of a 35-keV detector or a 3.9-MeV counter, or in terms of the marked change in character usually observed at the outer edge, 
or by the requirement that the flux change by more than a decade in a latitude interval of less than 0.5°. These four boundaries 
are plotted in Figure 14 as a function of local magnetic time. Physical interpretations have been given for the various boundaries. 
For example, comparisons with the auroral oval as given by Feldstein (1969) have shown that the boundary of the auroral oval 
coincides rather well with the so-called smooth boundary of Burrows and McDiarmid. 

A more sophisticated boundary definition has been used by Fritz (1968) in his analysis of Injun 3 data. His boundary is based 
on the degree of isotropy in the fluxes of electrons above 40 keV as measured with two detectors, one parallel to and one perpen- 
dicular to the local magnetic field line. The boundary between the high-latitude occurrence of isotropy and the lower latitude re- 
gion of anisotropy was usually sharply defined at local night but was rather broad or nonexistent at local afternoon during days of 
little geomagnetic activity. Fritz’s isotropic boundary was found to be in rather good agreement with the smooth boundary of 
McDiarmid and Burrows. The measurements of both McDiarmid and Burrows and of Fritz have indicated a narrow spike in the 
electron flux associated with the sharp boundary on the nightside. As the measurements of Imhof et al. (1969) indicate, the en- 
ergy spectrum associated with this spike may be much softer than that of electrons at slightly lower latitudes. Burrows and 
McDiarmid point out that the spike seems to represent an unstable region inside the nighttime sharp boundary. 


POLAR-CAP REGIONS 

Satellite-borne observations of solar protons at low altitudes over the earth’s polar caps have been reported by several groups 
of experimenters (Pieper et al., 1962; Akasofu et al., 1963; Stone, 1964 and 1969; Krimigis and Van Allen, 1967; Bostrom et al., 



Figure 14— A comparison of different outer radiation zone boundaries, as measured 
from Alouette 2 satellite (Burrows and McDiarmid, 1968). 
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1967; Paulikas et ah, 1968; Blake et al., 1968; Katz and Roth well, 1968; and Williams and Bostrom, 1969). Examples of energy 
spectra and intensities (Paulikas et aL, 1968 and Reagan and Imhof, 1969b) for the precipitating protons are shown in Figure 15. 

For several years it has been known that the latitudes of magnetic cutoff for solar protons deviate significantly from the pre- 
dictions of Stormer for a dipole field. Pieper et al. (1962) noted that the cutoff geomagnetic latitudes were of the order of 61° 
during the prestorm period; during the main phase of the storm the protons were detectable down to latitudes as low as 53° or to 
an £ -value of 3.0. Even the higher of these cutoff latitudes is lower than that of the classical Stormer theory. Stone ( 1 964) ob- 
served L5-MeV protons over the polar caps during a geomagnetically quiet period. The cutoff latitudes were 65° and 67° on the 
night and day sides, respectively, compared with a theoretical Stormer cutoff of 76°. Paulikas et al. (1968) have similarly observed 
lower latitude cutoffs than those predicted by Stormer. As is shown in Figure 16, the nightside cutoffs for protons of various en- 
ergies in the 1 to I GQ-MeV interval were at lower latitudes than calculated for a dipole field and in rather good agreement with the 
calculations of Reid and Sauer (1967) for a magnetic field model with the permanent existence of a geomagnetic tail. 

Solar protons have been observed to display pronounced spatial/temporal variations in flux over the polar caps (Reid and 
Sauer, 1967; Williams and Bostrom, 1967 and Blake et ah, 1968). This latitude structure greatly complicates the description of a 
given solar-flare event. As pointed out by Blake et al., a polar-cap average may not be an easily interpretable quantity. 


The simultaneous observations of solar particles from a low-altitude satellite in one of the polar regions and from a spacecraft 
in interplanetary space have been reported in several instances. Simultaneous observations have involved Explorer 7 and Pioneer 5, 
Injun 1 and Explorer 12, Mariner 4 and Injun 4, and Injun 4 and Explorer 33. In the latter set of measurements performed on the 
satellites Injun 4 and Explorer 33, located inside and outside the magnetosphere, respectively, Krimigis et al. (1967) found that 
0.5-MeV solar protons have access to the polar caps with a delay of 0.5 hour or less, as shown in Figure 17. These results were in- 
terpreted in terms of an open magnetospheric model with essentially immediate access from interplanetary space to the polar caps 
of the earth. Williams and Bostrom (1969) have reported general similarities between the flux-time profiles from Explorer 33 and 
from the low-altitude satellite 1963 38C, but have noted that significant differences in their time histories did occur. They have 
interpreted their observations in terms of the low-energy pro- 


ton access to the polar caps being controlled by a diffusion 
process. North-south asymmetries in the access of protons 
have been reported. Perhaps the most pronounced is the event 
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Figure 15— Energy spectra of protons observed over the polar 
caps at times of selected solar events. 


Figure 16— Nightside cutoffs as a function of energy. Calcula- 
tions of Reid and Sauer and of the classical Stormer theory are 
also shown (Paulikas et al., 1968). 
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Figure 17— Counting rates obtained with similar detectors onboard the Explorer 33 and Injun 4 satellites. 

Each of the points for Injun 4 represents the average counting rate during a polarcap pass (Krimigis et al, 

1967). 

of November 2, 1967, reported by Stone (1969), in which the access delay to the northern cap was about 20 hours compared with 
less than one hour in the southern region. 

Hie interpretation of the access of solar protons to the polar caps as a function of time, latitude, and energy, in terms of the 
various magnetosphere models is still an open question. The subject has been introduced primarily to emphasize the importance of 
low-altitude measurements of solar events. Here the term low-altitude applies to a much broader altitude range than considered in 
the radiation belts. In particular, for an isotropic influx of solar protons as often observed over the upper hemisphere (Paulikas et 
al., 1968) the low-altitude region is distinguished primarily by its location within the magnetosphere rather than by the presence of 
the atmosphere. However, at very low satellite altitudes, the proton environment may also be sensitive to steep flux variations 
with altitude for nonisotropic pitch-angle distributions. Such angular distributions have often been observed at positions near the 
magnetic cutoff. For studies of quasi-trapping phenomena, measurements at very low altitudes can indeed provide unique and im- 
portant information. There is clearly a great need for more extensive low-altitude coverage with instrumentation possessing sub- 
stantially greater sensitivity and better temporal and spatial resolution. 

In addition to isotropic and trapped pitch-angle distributions, the observation of energetic protons on magnetic field lines 
from the earth has been reported by Katz and Rothwell (1968). These protons were observed over the northern polar cap at an 
altitude of about 4800 km but not over the South Pole where the satellite altitude was 500 km. They have in terp reted the surpris- 
ing existence of such protons in terms of very large electric fields (~0. 14 volts/meter) parallel to the magnetic field lines. More 
pitch-angle measurements are probably needed before such phenomena can be well understood. 

In addition to protons, alpha particles and electrons have been observed over the polar caps in conjunction with certain solar 
events. The lack of experimental attention to these particles can be attributed partly to the lower intensities at comparable ener- 
gies in the case of alpha particles and in the case of energetic electrons to the fact that they appear to be accompanied always by 
energetic protons, but not vice versa (Vampola, 1969). According to Vampola, electrons of energy above 200 keV do not appear 
to be a common feature of the polar-cap region. However, at the time of certain solar events, electrons may appear with a reason- 
ably hard spectrum, as is shown in Figure 18. In addition to the spectral measurements of Vampola, the figure also includes some 
recent observations on April 13, 1969 by the Lockheed group (Imhof and Reagan, 1969b) from the polar-orbiting satellite OV1 18. 

From the limited observations to date the minority-group particles appear to behave differently than protons in several re- 
spects. Solar alpha particles in the polar-cap regions have given no evidence of durable trapping (Blake et al., 1968). In contrast to 
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the protons, energetic solar electrons have displayed no spatial 
structure over the entire polar zone (Vampola, 1969); that is, 
electrons of solar origin seem to have uniform access to the 
polar caps. 

CONCLUSION 

Clearly, many of the types of measurements performed to 
date should be carried out on a more regular long-term basis. 
For example, the inner-belt proton population should be fol- 
lowed in detail over a complete solar cycle. However, with a 
few notable exceptions the bulk of the low-altitude measure- 
ments have been performed on rocket or short-term satellite 
missions. Much information on the long-term behavior of 
trapped protons has been obtained through the use of nuclear 
emulsions. However, long-term counter data are clearly needed 
in addition, since emulsions are integrating devices and, there- 
fore, cannot generally provide good spatial resolution. As the 
Starfish population decays, it is of great importance to con- 
tinue monitoring the trapped electrons, particularly at the 
lower edge of the inner belt. Since the outer radiation belt en- 
vironment is very dynamic, there is a clear need for more com- 
plete simultaneous coverage from several satellites at different 
positions. Of course, there still remains an important role for 
vehicles other than satellites. For example, a rocket can pro- 
vide nearly simultaneous coverage over a broad range of L-shells at their magnetic equators or over a broad interval of altitude on a 
given magnetic field line. 

Several characteristics of the lower edge of the inner radiation belt (the short-trapping lifetimes, the low fluxes of trapped 
electrons, and the relatively small volume of space covered) have suggested that it may be an ideal region for increasing our knowl- 
edge of loss processes through the controlled injection of electrons by artificial means (Imhof et al., 1967). Such methods might 
include the placement of an accelerator (Christofilos, 1959 and Singer, 1962), reactor (Carpenter and Cohen, 1965), or radioactive 
source (Christofilos, 1 966), onboard a satellite or rocket. In contrast to the injection of electrons by the detonation of a nuclear 
weapon, these techniques provide for greater control of the energy spread and positions of injection. It has been calculated that 
early time fluxes of about 10 2 to 10 3 electrons/cm 2 -sec would be achieved during one pass of a low-altitude satellite at L « 1.15 
with an accelerator having a beam energy and current of 1.35 MeV and 100 j ua, respectively. Although these intensities are rela- 
tively low, they appear to be higher than naturally existing fluxes of electrons trapped at L » 1.15 and should provide adequate 
counting rates in a properly designed detector. By injecting a relatively small number of electrons at positions where the lifetimes 
are quite short, controlled trapping experiments can be performed without causing any long-term contaminations of the belts. 

In addition to the needs for better coverage of the low-altitude radiation environment with the type of instrumentation previ- 
ously employed, it is important to perform measurements with instruments of much greater sensitivity and higher spatial/temporal 
resolution. At low altitudes the fluxes are often quite small and the environment frequently changes very rapidly with position. 
Complete worldwide coverage with the required high-bandwidth instruments onboard the low-altitude satellites can now be 
achieved with improved data-handling capabilities. 



Figure 18—Energy spectra of electrons observed over the polar 
caps at times of selected solar events. 
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ABSTRACT 


Data from electron and proton detectors aboard satellite 1963 38C have been studied for the time 
period from September 1963 through December 1968. Detector energy ranges, in MeV, are E Q > 0.28, 
E e > 1.2, 1.2 < E p < 2.2, 2.2 < E p < 8.2, 8.2 25, and 25<E p <\ 00. The region to be dis- 

cussed is 1.2 < L < 3.0i? e over the Grange accessible to the satellite in an 1 100km circular, polar or- 
bit. For high-energy electrons (>1.2 MeV) the major feature has been the steady decay of the artificial 
belt. By the end of 1968, Starfish electrons had essentially disappeared and there was no indication of 
a natural source for electrons with energy >1.2 MeV. The low-energy (>0.28 MeV) electron flux has 
been near enough to the natural levels since early 1966 to reveal the effect of natural events; these 
electrons undergo significant changes throughout the inner zone in response to major magnetic storms. 
In particular, the events in September 1966, May 1967, and November 1968 produced large changes in 
the low-energy electrons and also in low-energy protons, although the protons responded over a more 
limited region of L- space. 


INTRODUCTION 

The temporal variations in the natural charged-particle population of the inner radiation zone have been studied by many 
authors (Pizzella et al., 1962; Pizzella, 1963; Mcllwain, 1965; Zmuda et ah, 1965; Imhoff et ah, 1967; Heckman and Nakano, 
1969; Bostrom and Armstrong, 1966; Filz and Holeman, 1965; Filiius, 1966; and Krimigis and Bums, 1969). The major perturba- 
tion to the inner zone produced by the Starfish high-altitude nuclear explosion has also been extensively studied (Van Allen, 1966; 
Hess, 1963; Mcllwain, 1966; Brown, 1966; Beall et al., 1967; Paulikas et ah, 1967; Pfitzer, 1968; and Bostrom et ah, 1968). A 
general summary of these observations would include the following statements. The inner radiation zone is a relatively stable re- 
gion in which significant changes in the high-energy particle fluxes are rare. Energetic protons are particularly stable; the perturba- 
tions which do occur are usually small and are limited to the outer region of the inner zone. High-energy (>1 MeV) electrons have 
a mean lifetime (1/e) of ^1 year in the inner belt but, apparently, are not a major constituent of the natural population. Low- 
energy electrons (<1 MeV) undergo large variations in intensity throughout the inner zone in response to major geomagnetic 
disturbances. 

This paper presents data on inner-zone trapped electrons and protons over a 63-month period from October 1963 through 
December 1968, from satellite 1963 38C. 

Satellite 1963 38C was launched into a circular, polar orbit at 1 100-km altitude on September 28, 1963. The satellite is mag- 
netically stabilized and contains an electron spectrometer and a proton spectrometer which measure locally mirroring particles. 

The instruments have been described in detail elsewhere (Williams and Smith, 1965; Bostrom et ah, 1967; and Beall et al., 1967) 
and pertinent characteristics of selected channels are summarized in Table 1. 

The coverage in B-, L - space obtained in November 1963 is shown in Figure 1 and is representative of the coverage during the 
lifetime of the satellite. The shape of the region is determined by the satellite orbit and the distribution of active ground stations, 
since there is no data-storage capability on board. The coverage, of course, goes on to very high L *s. All of the satellite data over 
its full lifetime have been scanned and 5-point (~5 sec) averages computed at 32 selected values of L , ranging from L- 1.2 to 
L~ 20. Nineteen detector counting rates, time, L, and B comprise each selected data point, and certain temperatures and voltages 
are retained for each pass. The passes through 1968 have been time ordered and the information presented in this paper is ob- 
tained from the selected, time-ordered data. The bars in Figure 1 show the B - ranges used in the time-history graphs to be shown 
later. The ^-restriction is necessary to reduce the scatter in the data. The i?-ranges associated with each L -value discussed herein 
are given in Table 2. 

An example of a computer-generated time plot is shown in Figure 2 (Beall, 1969). These are data for L = 1.25 and B between 
0. 16 and 0. 1 7 gauss, during 1964. Electron spectrometer counting rates are plotted with energy thresholds as shown. At this L 
measurements are made near the magnetic equator and the data are tightly grouped. Notice that there are two periods when the 
scatter in the data increases, Day 170 to 210 and from about Day 300 on. These are the times in 1964 when the satellite orbit was 
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Table 1 

Detector Characteristics. 


Channel 

Energy sensitivity (MeV) 

Geometric 

factor 

(cm 2 -ster) 

Electrons 

Protons 

Cl 

>0.28 

1.93 to 2.28; 
>153 

2.8 X 10~ 3 * 

C2 

>2.20 

14.4 to 14.5; 
>154 

2.8 X 10" 3 * 

(P1-P5) 

>0.7 

(low efficiency) 

1.2 to 2.2 

3.0 X 10~ 2 

(P2-P6) 

nil 

2.2 to 8.2 

3.0 X 10' 2 

P3 

nil 

8.2 to 25 

1.4 X 10" 2 

P4 

nil 

25 to 100 

1.4 X 10“ 2 



Table 2 

Range of 5-values as a Function of L. 


L -value 

Associated 5-range 

1.20 

0. 180 <5 < 0.190 

1.21 

0.175 <5 < 0.185 

1.22 

0.175 <5 < 0.185 

1.23 

0.170 <5 < 0.180 

1.25 . 

0.160 <5 < 0.170 

1.27 

0.160 <5 < 0.170 

1.30 

0.165 <5 < 0.175 

1.35 

0.170 < 5 < 0.180 

1.40 

0.180 <5 < 0. 190 

1.50 

0.190 <5 <0.200 

1.60 

0.195 <5 <0.205 

1.70 

0.200 <5 < 0.210 

1.80 

0.200 <5 <0.210 

1.90 

0.205 <5 <0.225 

2.00 

0.205 <5 <0.225 

2.20 

0.210 <5 <0.230 

2.40 

0.210 <5 <0.230 

2.60 

0.215 <5 <0.235 

2.80 

0.220 <5 <0.240 

3.00 

0.225 <5 <0.245 


in or near full sunlight, and the apparent reason for the scatter 
is an increase in the amplitude of the oscillation of the satellite 
about the magnetic field. The increased amplitude has also 
been observed by the onboard magnetometers. Notice that the 
scatter is only toward lower counting rates, as would be ex- 
pected if the detectors were measuring a wide range of local 
pitch angles instead of a narrow range near 90% In the follow- 
ing section ten-day averages of data such as these are shown. 

In all cases, the data were edited to remove these low points as 
well as some obviously noisy data. 

LONG-TERM OBSERVATIONS 


In Figure 3 the 5-year and 3-month time histories of elec- 
trons with energy >1.2 MeV for L-shells 1.2, 1.3, 1.4, 1.5, and 
1.6 are shown. These are 10-day averages of the counting 
rates. Bad points have been removed and dead-time correc- 
tions made. The main feature in the behavior of these particles 
is a monotonic decrease with time to a background counting 
rate that is due to energetic protons. The proton background 
level at each 5-value, derived from observations of the higher 
energy electron channels and based upon the procedure given by Beall et al. (1967) is shown. These levels differ only slightly from 
the values, given in the earlier paper, which were determined using only data through 1965. Of special interest is the fact that in 
no case is there any obvious response to the major magnetic disturbances which occurred during the satellite lifetime. This is in 
sharp contrast to the behavior of the lower energy electrons and implies that there is no strong natural source for 1.2-MeV electrons 


Figure 1 — 1963 38C { 8,L ! coverage in November 1963. This 
is representative of the coverage throughout the period from 
October 1963 through December 1968. Coverage extends to 

high L to larger B. 


*The procedure for obtaining flux from counting rate is given in Beall et al. (1967) and Bostrom (1969). The efficiencies for Cl and C2 are 
spectrum dependent but approximate values are ej - 0.77 and ~ 0.20 to 0.34. 
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Figure 2— Electron spectrometer counting rates versus time for L ~ 1.25 R e , 0.16 < B < 

0.17 gauss in 1964. 
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Figure 3— Ten-day average counting rates for C 2 {E e > 1.2 MeV) from October 1963 
through December 1968 at L - 1.2, 1.3, 1.4, 1.5, and 1.6 R e . The D gT scale indicates 
times of major magnetic storms. 

in the inner zone. The effective intensity threshold of our detector ranges from 10 4 to 10 5 electrons/cm 2 -sec-ster, depending on 
proton background and the electron spectrum. With that qualification, we would conclude that all electrons above about 1 MeV 
in the inner zone were injected by the Starfish event in July 1962. Pfitzer (1968) has observed inner-zone variations only for 
electrons below 690 keV. 

Curves of counting rate versus time, drawn through 10-day averages of the data from our lowest energy electron channel 
(threshold at 280 keV) are given in Figure 4. Comparison of this data at L- 1.2, 1.3, and 1.5 with that of Figure 3 shows that, 
whereas the high-energy electrons showed almost no structure, the low-energy electrons are strongly influenced by major magnetic 
storms throughout the inner zone. The Starfish component at this energy was undoubtedly large at L = 1.2 and 1.3 in 1963, but 
the intensity, even as low as L - 1.5, may well have been influenced by the September 1963 storm. This is clearly the case at 
L = 1.8 and 2.2. Storm-time changes show that the inner-zone electrons above 280 keV have been close to their (fluctuating) 
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Figure 4— Curves derived from ten-day average counting rates for C 1 (E e > 0.28 MeV) 
from October 1963 through December 1968 at L = 1.2, 1.3, 1.5, 1.8, and 2.2 R e . The 
D SJ scale indicates times of major magnetic storms. 


natural level of intensity since the beginning of 1966. In fact, at Z = 1.5, 1.8, and 2.2 a minimum in the intensity occurred in 
1966. The minimum for L - 1.2 and 1.3 occurred in 1967 but this may well have been because of the Starfish population at these 
/rvalues. Note that all data presented here are for /-values such that > 600 km and solar-cycle variations in atmospheric den- 
sity should not greatly affect the trapped radiation (Blanchard and Hess, 1964). The minimum is rather broad and appears to lag 
behind the minimum in the sunspot cycle by perhaps one to two years. Both the time lag and its /-dependence could be due to 
the presence of Starfish electrons. During this time period, major magnetic disturbances occurred in September 1966, May 1967, 
and November 1968 and, in each case, the low-energy electron population increased throughout the entire inner radiation zone. 
Smaller disturbances produce significant changes in the outer region of the inner zone; they probably produce small changes 
through the inner zone. This is examined in more detail in a later section. 

In the next four figures data on protons in four different energy ranges are presented. 

The 10-day average counting rates for protons with energy between 25 and 100 MeV at Z’s from 1.2- to 2.2-R e are shown in 
Figure 5. The maximum intensity in our orbit occurs at / ^ 1.27 and the flux (/j_) during the 1965-1968 period is 8.6 X 10 3 per 
cm 2 -sec-ster, compared with the value of about 2 X 10 3 from the AP 1 map produced by Vette (1966). Time variations are almost 
negligible. The September 23, 1963 storm produced a decrease in energetic protons at Z ^ 1.8, as reported by Mcllwain (1965). 
1963 38C was launched just after this event and the recovery is observed at Z = 1.8 and Z = 2.2. The only other major perturba- 
tion for these protons over this time period was in May 1967 and resulted in a decrease at / > 2. At Z = 2.2 the flux decreased by 
a factor of about three and one-half. Note that for a period of 4 years at L - 1.27 and 5 years at L- 1.2, 1.5, and 1.8, the count- 
ing rates have remained essentially constant with no obvious long-term trend associated with this period of the solar cycle. Data at 
Z = 1.27 through 1964 and at Z = 1.20 and Z = 1.50 through 1963 have been found to be contaminated by the intense Starfish 
electron fluxes and are, therefore, omitted. 

The data in Figure 6 are for protons between 8.2 and 25 MeV. Again five /-values beginning with 1.2 and extending to 2.6 
are shown. The others were chosen to show different regions while minimizing overlap. At this energy, there are very few rapid, 
large changes, although the September 1963 storm apparently had an effect and the May 1967 storm produced a decrease above 
1 = 2.2 which is clear at / = 2.6. However, in the 8- to 25-MeV range there is a slow but steady decrease over the 4 years, at / 
between 1.35 and 2.2. The slope of the line through the points at / = 1.35 yields a mean life of 7. 1 ± 1.4 years. At / = 1.2 there is 
essentially no change over this time period. 

The next lower energy, 2.2 to 8.2 MeV, where one begins to see considerably more activity is shown in Figure 7. Data for / 
between 1.5 and 3.0 are shown. The high-energy proton background has been subtracted, and this channel is insensitive to elec- 
trons. The major disturbances, when |D ST | exceeded 100 gamma, are shown across the bottom of the figure. There is a long-term 
decrease at all /-values shown and the effect of the September 1963 storm is not known. Except at / = 3, there are no large 
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Figure 5— Ten-day average counting rates for P4 (25 <£ p < 100 MeV) from October 
1 963 through December 1 968 at L = 1 .20, 1 .27, 1 .50, 1 .80, and 2.20 R g . 
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Figure 6— Curves derived from ten-day average counting rates for P3 (8.2 <£: <25 
MeV) from October 1963 through December 1968 at L = 1.20, 1.35, 1.90, 2.20, and 
2.60 R e . Representative ten-day average points are shown for L - 1 .35. 


changes associated with magnetic disturbances until May 1967. At L = 3 the response to major events may be a decrease as in 
April 1 965 and May 1 967, or an increase as in September 1966. The magnitude of the change is a factor of two to four and it 
takes place over a period of 20 to 40 days. The May 26, 1967 event obviously stands out as the largest magnetic disturbance in 
these 5 years, when D ST reached -370 gamma. It had the greatest effect on the inner-zone trapped particles during the period of 
observation. For these protons the largest enhancement occurred at L = 2.2 (a factor of nine) and this was accompanied by a de- 
pletion at L - 1.5 and at L - 2.6 by a factor of two. At L = 1.8 there was no change. These changes are examined further in the 
next section. 

In Figure 8, 1.2- to 2.2-MeV protons at L > 1.8 are shown. The vertical scale is expanded, and the high-energy proton back- 
ground has been subtracted. Although this channel is sensitive to electrons above 800 keV, the simultaneous data from the elec- 
tron spectrometer show that background from this source can contribute only at L ^ 2.6 and only during the first few days of an 
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Figure 7— Curves derived from ten-day average counting rates for (P2-P6) (2.2 <£ p < 
8.2 MeV) from October 1963 through December 1968 at L = 1.5, 1.8, 2.2, 2.6, and 
3.0 R e . The D sx scale indicates times of major magnetic storms. 
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Figure 8—Curves derived from ten-day average counting rates for (P1-P5) (1. 2<E„< 
2.2 MeV) from October 1963 through December 1968 at L- 1.8, 2.0, 2.2, 2.6, and 
3.0 R e . The D SJ scale indicates times of major magnetic disturbances. 


event. That is, the apparent lifetimes of energetic electrons are much shorter than those of protons and the temporal behavior of 
energetic electrons is quite different. We conclude, therefore, that at L- 1.8, 2.0, and 2.2, electrons do not contribute and that at 
L = 2.6 and 3.0 the electron contamination is small and of short duration (< 1 0 days). Here again a variety of responses to mag- 
netic disturbances are seen. In April 1965, there is a 50% increase at L = 2.6 and a decrease at L = 3 to about one-half the pre- 
storm value. Hie decrease at L = 3 is rather slow, taking about 75 days, with a further response to the storm on June 15th. In 
September 1966 the flux increased by a factor of three and one-half at L = 3, and by about 30% at L = 2.6. The response at 
Z < 2.2 was small or nonexistent until the May 1967 event, when the flux at Z = 2.2 increased by a factor of about 17. Z = 2.6 
showed no immediate effect and Z = 3, a rapid increase. At Z = 2 a gradual increase by about a factor of four occurred over a 
period of 1 00 days. 
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STORM-TIME VARIATIONS 


The primary purpose of this paper is the presentation of long-term observations on the inner-zone trapped particles. Each of 
the events associated with the major magnetic disturbances requires detailed study and correlation with other observations. How- 
ever, some aspects of the storm-time changes have been examined in slightly more detail than is shown on the 10-clay average 
graphs. These observations are presented in this section. 

The May 1967 event was, of course, the largest and produced major changes in the E e > 0.28 MeV electron population 
throughout the inner zone. A sequence of profiles in L of the Cl counting rate is given in Figure 9. 

Tire dashed line shows the situation prior to the September 1966 storm. The Day 130, 1967 profile shows that the long-term 
effect of the September 1966 storm was to enhance the region between L - 1.4 and 2.2. Day 140 is essentially identical with 130 
at L below 2.2, but the recovery from an earlier storm on May 3rd is seen at higher L. On Day 150, four days after the D gx maxi- 
mum, we find a large enhancement at higher L , 3 decades at L = 3, and also at low 5, an increase by a factor of five at L “ 1.2. 
However, a XL - 1.6, 1.7, and 1.8 the changes were relatively small. By Day 160, the region between L= 1.5 and 2.0 has been in- 
creased, and the region above L = 2 has begun to decrease. Finally, 40 days later on Day 200, L > 2.4 had returned to prestorm 
values; the flux at L~ 1.2 has returned to about twice its prestorm value; and the region between L = 1.5 and 2.4 is still well above 
prestorm levels by a factor of two to eight. At L = 1.5 the increase at lower L appears to be continuing, although the profile still 
shows a dip at L = 1.5. 

To examine these changes further, look at the storm-period data in more detail. Figure 10 was derived by means of all of the 
available data, which have been corrected for the 5-dependence. The 5-values to which the data were normalized are for L ~ 1.25 
and 1.3, B - 0. 18; for L = 1.2, 1.4, and 1.5, 5 = 0.19; for L = 1.6 and 1.7, 5 = 0.20; and for L - 1.8, 1.9, and 2.0, 5 = 0.21 gauss. 
There are two distinct regions of the inner zone in which particle injection and/or acceleration occur: (1) deep in the inner zone 
(1.2 1.4) where immediate and rapid increases are observed, and (2) the region of high L (L^ 1.9) where less rapid but 

larger percentage changes are found. In between is a transition region (1.5 *£L ^ 1.8) in which the increases are smaller and/or 
occur later. Although we are unable to determine very precisely the rise times, it is clear from the data that the increase at L = 1.2 
to 1.4 occurs in a few hours. This deep inner-zone effect is considerably smaller at L = 1.5 and 1.6, while at L = 1.7 it is not seen 
at all. Meanwhile, the much larger enhancement above L - 2 appears to move inward gradually. By Day 151 the intensity at 
L =]. 6 has begun to increase but there is no sign of this inward diffusion affecting L = 1.5. It is possible that the rapid increase at 
low L is the result of local acceleration since the detector is a threshold device. The high-L behavior is similar to earlier observa- 
tions of storm-time effects (Williams et al., 1968) wherein the 
largest and most rapid increase occurs at a particular value of L 
depending on the severity of the magnetic disturbance (in this 
case L » 2.4) and then the particle disturbance propagates 
both inward and outward in L. It is interesting, for instance, 
that the May 26, 1967 storm, which had such a marked effect 
on the inner zone, was not a particular!}/ noteworthy event in 
the outer zone. At L > 3.5 its effects were less notable than 
those of a number of much smaller storms. 

Another example of the rapid storm-time increases at low 
L was seen in November 1968. The clearly recognizable re- 
sponses at L = 1.2 to L - 1 .3 of electrons (E & > 0.28 MeV) to 
geomagnetic activity, as measured by XK p , are given in Figure 
1 1 . Other responses of the electron population might have 
occurred during periods of increased XK but due to sparse 
data they are less obvious than those shown. In fact, there 
seems to be a very close correspondence between XK p and the 
intensity of these electrons, but this topic will require more 
detailed treatment. Of primary interest here is the- electron 
response to the magnetic disturbance that began on Day 304 
and lasted through Day 309. Generally, both the rise and de- 
cay are slower at higher L. The uncertainty in the initial time 
of rise and the time of peak enhancement, shown in Figure 1 1, 
is ±1 day. The magnitude of the increase has a maximum at 



L (earth radii) 


Figure 9— Cl counting rate versus L profiles for Day 190, 1966 
(prior to September 1966 event), and for Days 130, 140, 150, 
160, and 200, 1967. The May 1967 storm occurred on Days 
145 to 146. 
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Figure 10— Cl counting rate versus time for 1.20 < L <2.0 R e in May 1967. D SJ 
values are plotted below the data. 


L = 1,21. Our detector does not permit us to determine whether we might be observing an energy-selective redistribution of the 
type reported by Imhoff et al. (1966) or some other acceleration mechanism. 

Hie May 1967 storm produced some major changes in the proton distribution as well. Hie magnitude of the changes is much 
smaller for protons than for electrons and these changes occur only in a limited region of space. The net effect of the storm is 
shown in Figure 12 where the prestorm and poststorm flux distributions for the four proton-energy ranges are presented. The flux 
from the AP 5 map produced by King (1967) is given for the 1.2- to 2.2-MeV range for comparison. The protons above 8.2 MeV 
were depleted at all U s shown, the depletion factor being larger for higher energy and higher L . The spectrum of the protons in- 
troduced at L = 2.2 apparently included some protons above 8.2 MeV, accounting for the small (15%) net loss of 8.2- to 25-MeV 
protons. For 2.2 to 8.2 MeV, there was little change at L < 2, an increase by a factor of nine at L ~ 2.2, an increase of less than a 
factor of two at L ™ 2.4, and a decrease at L > 2.5. For the lowest energy interval, the increase occurs over a wider range of L , 
but it still peaked at L = 2.2 and the only apparent decrease was at L = 2.6. If one examines a sequence of these cuts in time, the 
tendency is for the flux to decay at L > 2.2 and for the peak to broaden and increase at lower L. 


SUMMARY 

Data on trapped electrons and protons between L = 1.2 and 3.0 over a period of 5 years and 3 months, from September 1963 
through December 1 968 are presented in this paper. These data allow us to make the following statements. 

Protons above 25 MeV are extremely stable at L < 1.8. Above L = 2 they were depleted during the September 1963 and the 

May 1967 events. 

Protons between 8.2 and 25 MeV also respond only to exceptional magnetic disturbances but show a long-term decrease dur- 
ing this period. The apparent mean lifetime at L = 1.35 is 7.1 ± 1.4 years. 


32 




Figure 11— Cl counting rate versus time for 1.20<£<1.30 R e showing the low-L 
response to the November 1968 storm. K daily sums are plotted below the data. 


Low-energy protons (1.2 to 8.2 MeV) at L = 3 may be enhanced or depleted following a magnetic disturbance. The enhance- 
ment following the May 1967 storm was sharply peaked at L = 2.2 and accompanied by depletion at both higher and lower 
Z-values. The proton behavior strongly suggests a resonance phenomenon because the response is restricted in both energy and L, 

With regard to electrons above 1.2 MeV, in 5 years and 3 months we have seen no indication of response to magnetic pertur- 
bations at L < 2, there being only a monotonic decrease to a high-energy proton background. We conclude, therefore, that there 
is no source for 1.2-MeV electrons which can produce a quiescent or transient flux of more than 1G 4 to 10 5 electron$/cm 2 -sec~ster. 
The Starfish component at the end of 1968 was not measurable with our detectors. 

Electrons above 280 keV have been at or near their natural level of intensity since the beginning of 1966 for L > 1.3. There 
is some evidence that Starfish electrons persisted into early 1967 for LK 1.3. The response of electrons of this energy to a mag- 
netic disturbance is always an increased intensity. There are two distinct regions of response defined by immediate and rapid in- 
creases in electron intensity. The outer-region response is consistent with the previously reported outer-zone electron response to 
geomagnetic activity except that it is displaced into the outer edge of the inner zone for the very large magnetic disturbances that 
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Figure 12—Proton flux versus L for several energies showing the effect of the May 1967 storm on the 

proton distribution. 


were studied (peaking at L - 2.4 for the May 1967 event). The inner-region response is believed to be a previously unreported phe- 
nomenon. This response would have been masked in experiments prior to 1966 by the Starfish electrons. That it is a real phenom- 
enon is demonstrated by the close correspondence between electron intensity at L = 1.2 and 2K p (not presented in detail in this 
report) even for times of relatively low geomagnetic activity. These two regions are separated by a transition region, the inner edge 
of which is marked by an immediate but slow rise in intensity and the outer edge of which is marked by a delayed and slow rise in 
intensify. The inner edge of this transition region moves outward in L and the outer edge moves inward in L as the severity of the 
magnetic disturbance increases. As an example, in May 1967 we observed a rapid increase above L = 1.8 and below L = 1.6. The 
initial change was followed by an increase between L = 1.6 and L = 1.8 and subsequent decay at all U s. 


ACKNOWLEDGMENTS 

We wish to thank all of the many persons who have contributed to the long and fruitful life of satellite 1963 38C. In particu- 
lar, we gratefully acknowledge the assistance of Mr. C. J. Monahan, Miss B. A. Dodge, and Mrs. R. L. McQuiney in the analysis and 

presentation of the data. 

This work was supported in part by the National Aeronautics and Space Administration under DPR No. R2 1-009-004, and in 
part by the Naval Ordnance Systems Command, Department of the Navy, under contract NOw 62-0604-c. 


34 



REFERENCES 


Beall, D. S., “Graphs of Selected Data From Satellite 1963 38C for 1963-1967,” The Johns Hopkins Univ. Applied Physics Lab., 
Tech. Memorandum TG 1050 (5 volumes), April 1969. 

Beall, D. S., Bostrom, C. O., and Williams, D. J., “Structure and Decay of the Starfish Radiation Belt, October 1963 to December 
1965,”/. Geophys . Res. 72: 3403-3424, 1967. 

Blanchard, R. C, and Hess, W. N., “Solar Cycle Changes in Inner Zone Protons,” J. Geophys. Res. 69: 3927-3937, 1964. 

Bostrom, C. O., “Data User’s Guide to Satellite 1963 38C,” The Johns Hopkins Univ. Applied Physics Lab., Tech. Memorandum 
TG 1068, 1969. 

Bostrom, C. O., and Armstrong, J. C., “The Effects of the September 20-23, 1963 Magnetic Disturbances on the Inner Zone Pro- 
tons,” Trans. Am. Geophys. Union 47: 479, 1966. 

Bostrom, C. O., Beall, D. S., and Armstrong, J. C., “Trapped Particles in the Inner Zone, 1963 Through 1967,” (abstract) Interna- 
tional Symposium on the Physics of the Magnetosphere, Washington, D.C., September 3-13, 1968. 

Bostrom, C. O., Kohl, J. W., and Williams, D. J., “The February 5, 1965 Solar Proton Event. 1. Time Histoiy and Spectrums 
Observed at 1100 Kilometers,”/ Geophys. Res. 72: 4487-4495, 1967. 

Brown, W. L., “Observations of Transient Behavior of Electrons in the Artificial Radiation Belts,” in Radiation Trapped in the 
Earth’s Magnetic Field, ed. by B. M. McCormac, Dordrecht, Holland: D. Reidel Publishing Co., 1966, pp. 610-633. 

Fillius, R. W., “Trapped Protons in the Inner Radiation Belt,” J. Geophys. Res. 71: 97-123, 1966. 

Filz, R. C., and Holeman, E., “Time and Altitude Dependence 55-MeV Trapped Protons, August 1961 to June 1964,” / Geophys. 
Res. 70: 5807, 1965. 

Heckman, H. H., and Nakano, G. H., “Low- Altitude Trapped Protons During Solar Minimum Period, 1962-1966,” / Geophys . 

Res. 74: 3575-3590, 1969. 

Hess, W. N., editor, “Collected Papers on the Artificial Radiation Belt From the July 9, 1962 Nuclear Detonation,” / Geophys . 

Res. 68: 605-758, 1963. 

Imhoff, W. L., Cladis, J. B., and Smith, R. V., “Observations of an Energy-Selective Redistribution of Trapped Electrons in the 
Inner Radiation Belt,” Planet. Space Sci. 14: 569-577, 1966. 

Imhoff, W. L., Reagan, J. B., and Smith, R. V., “Long-Term Study of Electrons Trapped on Low L Shells,” / Geophys. Res. 72: 
2371-2377, 1967. 

King, J. H., “Models of the Trapped Radiation Environment Vol. IV: Low Energy Protons,” NASA SP-3024, 1967. 

Krimigis, S. M., and Bums, A. L., “The Time Dependence of Low Energy Protons and Alpha Particles in the Radiation Belts,” 
(abstract) International Symposium on the Physics of the Magnetosphere, Washington, D.C., September 3-13, 1968. 

Mcllwain, C. E., “Measurements of Trapped Electron Intensities Made by the Explorer XV Satellite,” in Radiation Trapped in the 
Earth’s Magnetic Field , ed. by B. M. McCormac, Dordrecht, Holland: D. Reidel Publishing Co., 1966, pp. 593-609. 

Paulikas, G. A., Blake, J. B., and Freden, S. C., “Inner Zone Electrons in 1964 and 1965,”/ Geophys. Res. 72: 201 1-2020, 1967. 

Pfitzer, K. A., “An Experimental Study of Electron Fluxes From 50 keV to 4 MeV in the Inner Radiation Belt,” Univ. of Minne- 
sota School of Physics and Astronomy Technical Report CR-123, August 1968. 

Pizzella, G., “On the Inner Van Allen Zone of the Earth,” Nuovo Cimento 19: 867-882, 1963. 

Pizzella, G., Mcllwain, C. E., and Van Allen, J. A., “Time Variations of Intensity in the Earth’s Radiation Zone, October 1959 
Through December 1960,”/ Geophys. Res. 61: 1235-1253, 1962. 

Van Allen, J. A., “Spatial Distribution and Time Decay of the Intensities of Geomagnetically Trapped Electrons From the High 

Altitude Nuclear Burst of July 1962,” in Radiation Trapped in the Earth’s Magnetic Field , ed. by B. M. McCormac, Dordrecht, 
Holland: D. Reidel Publishing Co., 1966, pp. 575-592. 

Vette, J. I., “Models of the Trapped Radiation Environment. Vol. I: Inner Zone Protons and Electrons,” NASA SP-3024, 1966. 

Williams, D. J., and Smith, A. M., “Daytime Trapped Electron Intensities at High Latitudes of 1 1 00 Kilometers,” / Geophys. Res. 
70: 541-556, 1965. 

Williams, D. J., Arens, J. F., and Lanzerotti, L. J., “Observations of Trapped Electrons at Low and High Altitudes,” / Geophys. 

Res. 73: 5673-5696, 1968. 

Zmuda, A. J., Pieper, G. F., and Bostrom, C. O., “Trapped Protons in the South Atlantic Magnetic Anomaly, July Through Decem- 
ber 1961,” 3. “Magnetic Storms and Solar Proton Events,” / Geophys . Res. 70: 2045-2056, 1965. 


35 




L0I6-TER! BEHAVIOR OF ENERGETIC INNER-BELT PROTONS 


H. H. Heckman 

Lawrence Radiation Laboratory, University of California 
Berkeley, California 

P. J. Lindstrom 

Space Sciences Laboratory, University of California 
Berkeley, California 

and 


G. H. Nakano 

Research Laboratories, Lockheed Missiles and Space Company 
Palo Alto, California 




ABSTRACT 


This report summarizes experimental results since late 1 962 of the temporal behavior of low- 
altitude trapped protons, E > 57 MeV. These data are interpreted in terms of the time-dependent 
continuity equation, using a source function that is empirically deduced from the data. The observed 
changes in the flux of 63-MeV protons between 1963 and mid-1969 are in agreement with those ex- 
pected from solar-cycle changes in the atmosphere. Preliminary calculations that include semiannual 
as well as solar-cycle atmospheric density variations indicate that beginning in 1967 semiannual changes 
in the 60-MeV proton flux should be observable at altitudes below about 350 km. Our flux measure- 
ments are consistent with this computational result. 


INTRODUCTION 

Since late 1962, we have monitored the low-altitude trapped protons, E > 57 MeV, using emulsion detectors recovered 
from oriented, polar-orbiting satellites. The span of nearly 7 years during which these observations have taken place constitutes a 
significant part of the current 1 1-year solar cycle. As part of this work, properties of the inner radiation belt protons have been 
examined in detail for a period concurrent with solar-minimum activity, from September 1962 through June 1 966, by Heckman 
and Nakano (1969). Throughout the solar-minimum period the characteristic feature of the energetic trapped protons was their 
high degree of stability. 

Coincident with increasing solar activity in mid- 1966, there is substantial evidence for a diminution in the proton flux at 
E - 63 MeV (Nakano and Heckman, 1 968). As of May 1969, the flux levels between 220- and 450-km altitude have decreased by 
more than a factor of two, relative to the solar-minimum period. It appears that the present solar cycle is very near the maximum. 
Consequently, the low-altitude proton flux is at, or near, its minimum value for this solar cycle. 

This paper reports briefly on the latest results on the solar-cycle variations in the proton flux. The altitude-flux profiles and 
differential energy spectra of trapped protons during solar minimum are reviewed. The second part of this report presents some 
preliminary computational results, using the time-dependent continuity equation (Singer, 1958 and Dragt, 1969), which is com- 
pared with experimental data. 


PROPERTIES OF TRAPPED PROTONS 
220- TO 450-km ALTITUDE 

At low satellite altitudes, the inner radiation belt is detected only in the South Atlantic anomaly. This region, centered at 
approximately 34° S by 34° W, is the site of anomalously low intensities of the geomagnetic field and, hence, high fluxes of 
trapped radiation. The geometry of the South Atlantic anomaly is given in Figure 1. The solid-line contours are Injun 3 isoflux 
contours (in particles/cm 2 -sec) for protons in the energy interval 40 < E < 1 10 MeV at an altitude of 400 km (Valerio, 1 964). The 
contours of constant L , computed by using the 48-coefficient Jensen and Cain geomagnetic field model (1962), are shown as 
dashed lines. The northward trajectories traversing the anomaly region are those of a typical 4-day satellite flight at 400-km alti- 
tude and 75° orbit inclination. Demonstrated here is the uniform sampling of the particle fluxes in the anomaly, as well as the 
limited range of the Z-parameter. The magnetic field intensity, B y varies from 0.21 to 0.25 gauss over the range of L -values shown. 

The inner-belt protons exhibited a high degree of stability during the recent solar-minimum period, September 1962 through 
June 1966 (Heckman and Nakano, 1969). Some of those temporally stable features of the inner-belt protons referred to in the 
second part of this report are illustrated in Figures 2 through 5. The altitude dependence of the omnidirectional flux of E = 63 
MeV protons is presented in Figure 2. We plotted the average daily flux, J (in protons/cm 2 -day-MeV), normalized to a 90° orbit 
inclination, observed for each satellite flight, versus the average minimum mirror-point altitude, S min . This flux refers only to those 
protons detected during northward passes of the satellite in the anomaly (see Figure 1). To remarkable precision, the altitude pro- 
file of the omnidirectional flux of 63-MeV protons can be fitted by a simple power-law function of the form /(65 MeV) ec h 1 !-. 
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Figure 1-The South Atlantic anomaly. The heavy-line con- 
tours are Injun 3 isoflux contours for protons, 40 < E < 1 10 
IVIeV, at 400-km altitude. Dashed lines are contours of L from 
the Jensen and Cain 48-coefficient field model. The northward 
traversals of a satellite, 75° orbit inclination, during a 4-day 
flight are shown as arrows. 



hmin (km) 


Figure 2— Omnidirectional flux of 63-MeV protons versus aver- 
age minimum-mirror-point altitude during solar minimum, 
November 1982 to June 1966. Experimental data are shown 
as open circles; the least-squares fit to these is denoted by J c . 
Computed Injun 3 daily fluxes and least-squares fit are shown 
by the solid points and dashed line. 


The experimental data are shown as open circles and are fitted 
to a least-squares power law, J c « h^ ±om , to an accuracy of 
±7.6% (SD). The computed daily fluxes for each flight, de- 
rived from the Injun 3 B- , L -proton-flux contours and the sat- 
ellite ephemerides, are shown as dark circles. The least-squares 
fit to these data yields an exponent n - 4.05 ± 0.04. Although 
the differences between the slopes of the observed and com- 
puted altitude-flux profiles are statistically significant, we be- 
lieve the Injun 3 data (recorded between December 24, 1962 
and September 28, 1963) reproduce well our solar-minimum 
data, clearly within the stated accuracy (20%) of the Injun 3 
data. 

Five measurements of the differential energy spectra for 
protons, E > 57 MeV, between June 1963 and June 1966, are 
given in Figure 3. The data are normalized to an altitude h min = 
375 km, assuming the omnidirectional flux varies as in- 
dependent of energy. The errors of measurement are consist- 
ent with the spread in the data points. There is, therefore, no 
evidence for a change in the spectral shape of the protons dur- 
ing the solar-minimum period. In f^act, the energy spectrum 
obtained by Heckman and Armstrong (1960) (shown as a 
dashed line, normalized at E = 63 MeV) also exhibits the same 
spectral shape. Thus, within the observational errors, the pro- 
ton energy spectrum has been invariant for almost a complete 
solar cycle. 

In what follows, the spectral shape of the spectrum be- 
tween 50 and 90 MeV is needed. For this energy interval, we 
take the spectrum to be dJ(E)/dE E~ 0J2 . This energy spec- 
trum was first derived by Freden and White (1960) and is in 
agreement with the observed spectrum for 50 < E < 90 MeV. 

The temporal behavior of the omnidirectional proton flux 
at 63 MeV, between September 1962 and May 1969, is pre- 
sented in Figure 4. Plotted as a function of time is the quantity 
J(t)/J c , the ratio of the 63-MeV omnidirectional proton flux at 
time, t, to the least-squares flux, observed during the solar- 
minimum period (Figure 2). The data are grouped into three 
altitude ranges: (a) 220 < /z min < 300 km, (b) 300 < h min < 
400 km, and (c) 400 < h min <455 km. Within the experimen- 
tal errors (the 5% statistical errors are indicated for the 300- to 
400-km interval) we conclude that the flux of 63-MeV protons 
was, for all practical purposes, in a steady state at all altitudes 
between 220 and 455 km during the period November 1962 to 
June 1966. 

Variations in the proton flux before and after the quies- 
cent period are readily apparent. The high value of the flux 
observed in September 1962, the first measurement, can be 
attributed to the effects on the inner-belt protons of the Star- 
fish nuclear detonation in July 1962 (Filz and Holeman, 1965). 
Near mid- 1966, the proton flux began to decline and as of 
May 1969, had decreased to about 0.40 of its solar minimum 
value. 
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CONTINUITY EQUATION: COMPUTATIONS 
Solar Cycle Variations 

It is these experimental results that we wish to discuss in 
terms of the time-dependent continuity equation (Singer, 1958 
and Freden and White, 1 960) 

dN _ d _ „ __ dE , . 

— = S-—vN( E )- - pavN , (1) 

where 

N = particle density in particles/cm 3 -MeV 
and 

S - source strength in particles/ cm 3 ~MeV~sec. 

The second and third terms on the right are, respectively, the 
rates of particle loss due to energy loss by ionization and nu- 
clear interactions in the residual atmosphere. For the atmos- 
pheric density averaged over the particles 5 trajectory, p, we 
used the results of Heckman and Brady (1 966). In our present 
calculations, the results of Heckman and Brady, which are 
based on the Harris and Priester (1962 and 1963) model at- 
mosphere, are reformulated in terms of F l0 7 (in !G~ 22 W/m 2 - 
cps), the 10.7 solar flux. Furthermore, we have augmented 
the atmospheric density calculation to permit the evaluation of 
the effective atmospheric densities as a function of F 10 7 as 
well as E , the proton energy. The latter correction enters be- 
cause the atmosphere traversed by a trapped particle is depend- 
ent on its cyclotron radius, hence energy, for a given minimum 
mirror-point altitude of the particle’s guiding center (Heckman 
and Brady, 1966). 

The experimental data cited above, when incorporated 
into solutions of Equation 1, should give some insight as to the 
source strength, altitude dependence, and the temporal behav- 
ior of the trapped protons. Explicitly, the data state that, in 
the energy interval 50 <E < 90 MeV, 

dN 

-—jj ~ 0 (between 1963.0 and 1966.0) , 

(250 < h min < 450 km) , 

and 

N(E) oc E~ 0J 2 . 

Given the time-dependent effective atmospheric density, 

^mm ’ -^10 7 ^’ w ^ ere ^ 107(0 is monthly averaged 
10.7-cm solar flux, our objective is to solve for the energy and 
altitude dependence of the source function, S(E, h), from 
which the particle density N(E) can be computed as a function 
of time and altitude. 

To prepare Equation 1 for integration, we follow Freden and White (1960) and express the velocity, v, and rate of energy 
loss, dE/dx, in terms of power-law relationships in energy, E, for the interval 50 <E < 90 MeV. Given a source strength, S, and a 
proton density at time t , N{E, t), then, at time t + At, 
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Figure 3— Differential energy spectrum for protons, £>57 
MeV, 1963 to 1966. Dashed curve is a smoothed fit to the 
spectrum measured in 1960 by Heckman and Armstrong. 
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Figure 4-Temporal behavior of the 63-MeV omnidirectional 
proton flux, September 1962 to May 1969. Plotted are the 
ratios J(t)(J c versus time, where J c is the least-squares fit to 
the data for the period November 1962 to June 1966. 
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N(E, t + At) = N(JE, t) + Af . 


( 2 ) 


Integration errors in Equation 2 are less than 0.1% for step size At of 1 to 2 hours. 

Is there an S(E,h) such that an approximate steady-state solution is given by Equation 1 between 1963.0 and 1966.0, inde- 
pendent of E and hi Y es, there is. As a particular example, S(h) obs , for 60-MeV protons, that yields flux values that are practically 
constant during the solar-minimum period, 250 < h min < 450 km, are given in Figure 5. The errors assigned to *S W obs represent 
the uncertainty in the source function that arises from the approximate scatter of ±10% in the flux measurements about their 

mean value for the period 1963.0 to 1966.0. 

The significant result is that, between 250- and 450-km altitude, the source function is not constant (as one might expect 
were the proton source solely due to albedo-neutron decay), but decreases with increasing altitude according to Y(&) obs a 
Also presented are the altitude dependences of TV (at 60 MeV) ex h^ n (observed) and p /f 7 * 0 (computed). As noted by Cornwall 
et al. (1965) the quantity S/Np is a constant under equilibrium conditions. Our results are in good agreement with this condition. 
The numerical values of S and TV need be relative only, because it is the ratio S/N that is important here. However, to be somewhat 
realistic, we have adopted for TV (at 60 MeV) at 1963.0 the proton density obtained from the Injun 3 data for h min = 350 km, 

L = 1 .40 earth radii. 


In their examination of the changes in 2?-, Y-space with time, Lindstrom and Heckman (1968) pointed out that the minimum 
mirror-point altitudes in the South Atlantic anomaly are decreasing at a rate of h^l km/year owing to the secular change of the 
earth’s magnetic field. Such a drift in the minimum mirror-point trajectories in the anomaly produces an apparent source of pro- 
tons of the amount S B = h(dN/dh ), where h is the inward drift 



ALTITUDE, h min (krn) 


velocity and dN/dh is the altitude gradient of the proton flux. 
This apparent source is also shown. We note that its magnitude 
does not significantly contribute as a proton source at the 
lower altitudes, comprising only about 2% of 6' obs at 250 km. 
However, because of its strong dependence on altitude (that is, 
S B varies as & 3 * 7 ) the S B source becomes quite significant for 
higher altitudes, and can account for nearly one-half of S abs 
for altitudes above 450 km. It should be stressed that the ap- 
parent source, S B , arises solely from the dynamic behavior of 
the geomagnetic field and is independent of any assumptions 
as to nonadiabatic behavior in the particle motion. 

The results of computations on the solar-cycle variations 
in the flux of 60-MeV protons are presented in Figures 6 and 
7. Figure 6 shows the temporal behavior at 250, 350, and 450 
km at L = 1.40 between January 1963 and May 1969. Plotted 
versus time are the computed ratios 1963.0). The exper- 
imental data (Figure 4) are superimposed on the computed 
curves for comparison. 

The good agreement between the solutions of the conti- 
nuity equation and our experimental data is most encouraging. 
Between 1963.0 and 1966.0, the empirical source function, 
Sobs’ re P r °duces well the observed constancy of the proton 
flux. As solar activity developed after 1966, the computations 


show the characteristic differences in the behavior of 60-MeV 


Figure 5—Aititude dependence of source function of 60-MeV 
protons, S obs , required to account for a constant flux level ob- 
served between 1963.0 and 1966.0. The altitude dependence 
of the measured proton density, N{h), and computed effective 
atmospheric density, p{h), for 1963.0 are also shown. The 
source S & is the apparent source of 60-MeV protons arising 
from the lowering of the mirror-point trajectories owing to the 
secular variation of the geomagnetic field. 


protons as the altitude and, hence, lifetimes, increase. Since 
mid- 1 967, the calculations indicate that the proton flux at 250 
km has sustained a nearly constant minimum value, coincident 
with the plateau in activity during solar maximum. That is in 
contrast to the fluxes at 350 and 450 km, which were still de- 
creasing in early 1969. Such an effect accounts for the appar- 
ently high values of J/J c observed in our last two flights (h min = 
224 and 299 km). 
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Figure 6— Computed ratios of J(t)U{ 1963.0) at 60 MeV and minimum-mirror point 
altitudes of 250, 350, and 450 km, L = 1.40, between January 1963 and May 1969, 
given source S Qbs . Flux data from Figure 4 are also shown. 


Starfish Redistribution 

Filz and Holeman (1965) have reported that an abrupt increase in the low-altitude proton flux followed the Starfish nuclear 
detonation on July 9, 1962. At 350 ± 50 km, their data indicate the proton flux was increased by about seven times over its pre- 
bomb level by this event, owing to a nonadiabatic redistribution of the low-altitude protons. On the basis of Filz and Holeman' ’s 
work (1965) we would conclude that the proton flux decayed within 1 year to a constant level, which was some two and one-half 
times the prebomb level. It is evident, therefore, that the temporal behavior of the inner belt, which we have been studying, has 
likewise been artificially perturbed. To examine this further, we have carried out the integration of Equation 1 for the flux 
/(O 60 MeV ^min = 350 km beginning 1953.0, under the assumption that S Qhs is time independent. Two initial arbitrary flux val- 
ues were assumed: (a) /( 1953.0) = 9.2 protons/cm 2 -sec-MeV, equal to the 60-MeV flux observed in 1964, one full solar cycle later, 
and (b) /(1 953.0) = 3.2 protons/cm 2 -sec-MeV. 

The results of this computation are presented in Figure 7. Solutions (a) and (b) converge to a unique solution after 1958.5, 
yielding a J(t) that is independent of initial conditions. The anticorrelation of the proton flux with the atmospheric density, the 
latter shown by the dashed line is well illustrated. The computed curves of J(t) from 1958 through 1968 display how flux changes 
at 60 MeV would have occurred had there been no Starfish detonation and, hence, proton redistribution. Our interpretation of 
the sequence of events portrayed is that, on July 9, 1962, the Starfish detonation caused a rapid increase in the proton flux, and, 
after an apparently short-lived transient period, settled down to a 3-year, steady-state condition. This steady-state phenomenon 
was fortuitously produced when the natural solar-cycle increase in the proton flux was just compensated by the decay of the Star- 
fish perturbation. The difference between the Starfish and no-Starfish calculations is, in fact, an exponential decay curve. The 
particle data denoted by circles are from Filz and Holeman (1965); the data of this experiment are shown as X’s. All the data are 
restricted to h min = 350 ± 50 km. The two sets of data points are independently normalized to the computed, steady -state condi- 
tion for 1963.0 to 1966.0. We believe it is significant that the relative changes in the ( 5 7-MeV) proton flux observed by Filz and 
Holeman at 350-km altitude, before and after the Starfish detonation, are in agreement with calculations that were, in turn, based 
on our independently measured post-Starflsh data. Finally, we note that the effects of Starfish have now virtually disappeared for 
altitudes of 350 km and less. Henceforth, the inner-proton belt will be (one hopes) in its natural state at these low altitudes. 
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Figure 7-Computed proton flux v/(t ) 60 MeV at/? min = 350 km, beginning 1953.0, assum- 
ing source S Qbs is constant, and initial flux values 7(1953) = 9.2 and 3.2 cm" 2 sec" 1 
MeV" 1 . Unique solution is obtained at this altitude after 1958.5. The Starfish redistri- 
bution on July 9, 1962, the 3-year steady-state, 1963 to 1966, and the abrupt solar- 
cycle decrease in the proton flux after mid-1966 are reproduced in these calculations. 

The data denoted by open circles are from Filz and Holeman; the x's, this experiment. 

The data points are limited to h m]n = 350 ± 50 km. Solar-cycle behavior of atmosphere 
at 350-km altitude is given by the dashed curve. 

Semiannual Variations 

Recent satellite atmospheric-drag results have revealed that during the course of our particle measurements, the semiannual 
and solar-cycle density variations in the atmosphere have been comparable in magnitude at altitudes between 191 and 1030 km 
(King-Nele et aL, 1968 and 1969 and Cook and Scott, 1967), The semiannual variation has yet to be explained, although the effect 
has been well documented. It is adequate for our purposes that, to the extent it has been studied over an 8-year period, the semi- 
annual effect is a highly stable, global phenomenon (Jacchia et al., 1969). Irrespective of altitude and local time, the atmosphere 
exhibits in-phase world-wide density variations, the minima occurring on about January 26 and July 25, and the maxima on April 
1 and November 1. Although the periods between minima and maxima are not equal, the dates on which they occur have been 
determined to within an accuracy of several days. The minimum-to-maximum density ratios affected by the semiannual variation 
are typically 1:2. The observed amplitudes of the density variations appear to be related to temperature variations in the thermo- 
sphere, the latter being correlated with the intensity of the 10.7-cm solar flux. 

It is quite obvious, then, that when proton lifetimes are less than about 6 months, observable semiannual changes in the flux 
of inner-belt protons should occur. To investigate this statement quantitatively, we have attempted to include in our atmosphere 
model a realistic description of the semiannual density variations. To do so, we have referred to the works of King-Hele et al. 

(1968 and 1969) and Cook et al. (1967), who have, since 1964, examined the semiannual variations at altitudes 191, 470, 480, 
and 1 130 km. We have collected these data and presented them in Figure 8. To reveal the semiannual variation, diurnal and solar- 
activity effects have been removed from the data. The residual semiannual variations were then reduced to an F l0 7 that was char- 
acteristic of the time during which the atmospheric density measurements were made. The dashed lines through the data are the 
analytical functions we use to represent the semiannual effect. 

We now have a description of a (bounce- and drift-averaged) Harris and Priester atmospheric model that includes both solar- 
cycle and empirical semiannual time variations. The flux changes for 60-MeV protons predicted by the continuity equation with 
this model atmosphere are given in Figure 9. Upon comparing this figure with Figure 6, we see that the semiannual density changes 
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Figure 8— Semiannual density variations in the atmosphere, reproduced from King-Hele 
et al. and Cook and Scott. Dashed curves are analytic functions used to represent the 
semiannual density changes in the computations. 



TIME (yr) 


Figure 9-Calculated J(t)/J( 1963.0) at 60 MeV at h . = 250, 350, and 450 km when 

min 

both solar-cycle and semiannual density changes are incorporated into the continuity 
equation. Again, the experimental data are shown for comparison with the theory. 
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are clearly present at 250 km, but not so at 350 and 450 km. There are two salient features of these computations that are in 

accord with the experimental data. 

First, between 1963.0 and 1966.0, the semiannual changes in the proton fluxes are small; that is, they are typically equal to 
or less than the experimental errors. The computed results continue to exhibit a very stable flux level during solar minimum. 

Second, as solar activity increases, the average density of the atmosphere increases and the proton lifetimes decrease. Conse- 
quently, the response of the trapped protons to changes in the atmosphere becomes more pronounced. This is exemplified by the 
increasing semiannual oscillations in the proton flux since 1967. The observed amplitudes of the fluctuations in the proton-flux 
measurements since 1967 are consistent with this interpretation. 

Because of the preliminary and sensitive nature of the calculations to our atmospheric and magnetic field models, we do not 
attribute any significance to the fact that the variations in the 350 ± 50-km data appear to be characteristic of an altitude slightly 

less than 350 km. 

Thus far, a number of instances in which the experimental data are accounted for by theory have been discussed. In contrast 
to this, an interesting, and perhaps quite important, observation can be made by closely examining Figure 9. It is that the observed 
flux changes are not in phase with the computed semiannual flux variations. Taken at face value, the flux changes appear to be out 
of phase with the solutions of the continuity equation by about 5 weeks. Surprisingly, this phase difference is such that the proton- 
flux variations are nearly in phase with the semiannual density changes in the atmosphere. The significance of this observation is 
certainly not clear at this time. We wish, however, to make the following observation. The principal unknown in the theory of the 
energetic, inner-belt protons is the source. We have treated the source S(E 9 h) 9 therefore, as an empirical quantity to be determined 
by experiment. Using this approach, we find that this empirical source, N obs , decreases with increasing altitude. The only other 
physical entity that decreases with altitude is the density of the atmosphere. The implication, then, is that there is an effective 
source of trapped protons related somehow to the atmosphere, indicative of pitch-angle scattering. 


SUMMARY 

On the basis of our experimental studies on the temporal behavior of low-altitude, inner-belt protons, and their interpretation 
in terms of solutions of the time-dependent continuity equation, using an empirically determined source function, we make the 

following conclusions. 

The observed temporal changes in the flux of 63-MeV protons since 1963, 220 < h min < 450 km, are in agreement with those 
expected from the solar-cycle changes in the atmosphere. 

The computed time-history of the proton flux at 60 MeV before and after the Starfish detonation is in excellent agreement 
with the observations of Filz and Holeman, as well as with our post-Starfish data. 

For altitudes of about 250 km, semiannual density changes in the atmosphere should be observable in flux changes of pro- 
tons, E ^ 60 MeV. That a phase difference may exist between the observed and computed semiannual changes in the proton flux 

is unexpected and needs to be explained. 

Pitch- angle diffusion, which may be an effective source strongly correlated with the atmospheric density, must now be in- 
cluded in the theory' in order to see whether or not it can account for the altitude dependence of the source function implicated 

by our experimental observations. 
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INTRODUCTION 


The synchronous orbit is now perhaps the region of space most densely populated by satellites. The obvious advantages ob- 
tained by communications and meteorological instrumentation placed in this orbit have resulted in a satellite population explosion 
at the synchronous altitude. More than 50 satellites had been placed into synchronous or nearly synchronous orbit as of late 


1969.* 


The requirements for engineering data about the synchronous radiation environment were realized at an early date and the 
now historic electron flux compilation by Vette and Lucero (1967) resulted. It is interesting to note that no data from satellites in 
a truly synchronous orbit were used in the Vette/Lucero compilation since no synchronous satellite was instrumented with particle 
radiation detectors prior to 1 966. Data obtained at or near the synchronous orbit from satellites in elliptical orbits were used to 
construct the model electron environment. In a subsequent study, King (1967) presented a model of the low-energy trapped pro- 
ton environment for the outer zone, including the synchronous orbit. 

The first satellite, in synchronous orbit, equipped with an array of radiation detectors was the ATS 1 launched in December 
1966. Subsequently, other satellites with complements of particle instrumentation have been placed into this orbit. A picture of 
the synchronous particle population and the changes that this population undergoes with changing magnetic activity can now be 
constructed, based on data from geostationary satellites. Table 1 lists the satellites which are now in orbit and have contributed 
(or will, hopefully, contribute) to the picture of the particle environment at synchronous altitude. The amount of available data 
regarding the particle environment at this altitude is now immense; the present report will consider, in the main, data obtained by 
the Aerospace Corporation omnidirectional spectrometer on ATS 1 . 

The synchronous orbit is also very interesting from the point of view of radiation belt physics. Some separation of spatial 
variation from temporal effects can be obtained because the satellite samples only a limited range of the magnetic field. Field lines 
connect the synchronous region with the auroral zone and experiments correlating measurements on the magnetic equator with 
the foot of the field line can be readily done (Parks et al., 1968). Solar protons penetrate readily to this altitude and the question 
of access of solar particles to the polar regions is intimately connected with solar particle access to the synchronous altitude. 

Because the synchronous orbit coincides approximately with the center of the ring current (Frank, 1968), the particle popu- 
lation exhibits extremely dynamic behavior and construction of an average picture of the environment requires an almost statistical 
approach. Such an approach was initially used by Vette and Lucero to order the energetic electron data and this approach, where 
possible, has been adopted in this paper. 

While most of the time a synchronous satellite is comfortably within the domain of closed field lines, the synchronous orbit 
approaches the trapping boundary on the night side. During magnetically active periods, plasma intruding from the tail of the 
magnetosphere may reach down to these altitudes (Freeman, 1968 and Freeman and Maguire, 1967). This causes perturbations of 
the field lines on the night side of the earth (the field may sometimes become almost radial) and, therefore, very significantly mod- 
ifies the trapped particle population. 

In several instances on record, changes in the solar wind momentum were sufficiently large to compress the boundary of the 
earth’s magnetic field to a radius <6.6 earth radii. At those times, of course, the energetic particles disappeared and were replaced 
by low-energy particles characteristic of the magnetosheath. 

Solar protons have surprisingly free access to the synchronous region even during times of magnetic calm (Lanzerotti, 1969). 
During magnetically active periods, solar protons well below 1 MeV can enter the synchronous region with essentially full effi- 
ciency. Thus, solar particles also form an important constituent of the particle environment. 

In the subsequent sections of this paper, the behavior of the population of trapped energetic electrons and the penetration of 
solar protons to the synchronous orbit will be discussed in detail. Although present only temporarily at the synchronous altitude 
following solar proton events, solar protons are an important element in the total radiation picture near the period of maximum 
solar activity (Lanzerotti, 1969). 

* TRW Space Log. 
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Table 1 

Scientific Satellites at Synchronous Altitude. 


Satellite 

Launch 

ATS 1 

December 1966 

OV2 5 

September 1968 

ATS 5 

August 1969 


The limitations of the subject matter of the present re- 
port (high-energy electrons and protons) should not obscure 
the underlying unity and interrelation of the behavior of the 
trapped particle population. It is particularly important to re- 
member that the bulk of the particle energy density resides in 
the low-energy (approximately tens of keV) population of 
trapped protons and electrons. The behavior of the trapped 
energetic electrons and solar protons merely reflects the inter- 
action of the low-energy ring-current particles and the geo- 
magnetic field; energetic electrons and solar protons do not 
contribute much energy to this particle-field interaction. 


THE MAGNETIC CONTAINER* 

Central to the discussions of the particle population and the spatial and temporal variations of these particles is an understand- 
ing of the configuration of the geomagnetic field which confines and guides the motion of the charged particles. Figure 1 is a qual- 
itative picture of the magnetosphere and shows the location of the synchronous orbit (at 6.6 R e ) in relation to the main features of 
the magnetosphere. As can be seen, the earth’s magnetic field at high altitudes loses cylindrical symmetry and dipolelike character 
because of the distortion of the terrestrial field caused by the solar wind. A satellite at the synchronous orbit sees a diurnal varia- 
tion (in the absence of temporal variation) in the strength and orientation of the magnetic field as well as in the trapped particle 
population. The loss of cylindrical symmetry has an important effect upon trapped particle motion in that the parameter L be- 
comes a function of pitch angle, so that the longitudinal drift path is now pitch-angle dependent. The families of drift paths for 
two extremes of equatorial pitch angle are shown in Figure 2 in relation to the synchronous orbit. A satellite at 6.6 R e samples a 
range of particle-drift shells which, for the particular magnetospheric parameters used in the calculation of the curves of Figure 2, 
have equatorial crossings from 6.1 to 7.6 R e . Measurements of omnidirectional flux such as form the bulk of this compilation, 
therefore, are representative of a range of radial distance (Roederer, 1967, 1968, and 1969). 

An important consequence of the asymmetry of the magnetosphere is that particles can be locally trapped yet be unable to 
execute a complete longitudinal drift period. Such particles are called quasi-trapped, and their mirror points are located in the 
quasi-trapping regions of Figure 1 . Particles mirroring at high latitudes near local noon are connected with open lines in the tail. 
Conversely, particles mirroring at low latitudes near the midnight meridian will drift out the flanks of the magnetosphere and fail 
to complete a longitudinal circuit. The particle angular distribution which results is illustrated (in a qualitative manner) in Figure 
3. A drift loss cone appears, in addition to the usual bounce loss cone associated with particle loss into the atmosphere. In general, 
the drift loss cone can be much larger than the small (about 3°) bounce loss cone and thus modify the angular distributions in an 
important way. Since the synchronous orbit often enters the region of quasi-trapping, the existence of the drift loss cones must be 
explicitly considered. White et al. (1969) have experimentally confirmed the theoretical predictions by detecting the existence of 
a proton drift loss cone. 

The model of Mead and Williams (1965) used to calculate Figure 2, has proved to be a good representation of the quiet-time 
magnetic field at the synchronous orbit. Roederer et al. (1968) using ATS 1 magnetometer data have devised a procedure to adjust 
the parameters contained in the Mead-Williams model and in this way analytically reproduce the observed field components and 
the diurnal variation in these field components. The quiet-time magnetic field can thus be reasonably well related to a field model. 
(This model neglects such refinements as the diurnal and seasonal variation of the aspect of the earth’s magnetic field axis with re- 
spect to the solar wind.) In effect, then, extended measurements of magnetic field components at the synchronous altitude amount 
to synoptic monitoring of the shape and configuration of the entire magnetosphere (Roederer, 1969). 


ENERGETIC ELECTRONS AT THE SYNCHRONOUS ALTITUDE 


Typical Data 

To set the background for subsequent discussions, we present in this section several sequences of data that illustrate the vari- 
ability of the electron fluxes at synchronous altitude. The data show the transition from magnetically quiet to disturbed conditions 

* An extensive review of the properties of the magnetic field at the synchronous altitude can be found in a report by Coleman (1968). 
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Figure 1 — Location of synchronous orbit (dot at 6.6 R e ) in a 
Mead-Will iams model of the magnetosphere. Particles mirror- 
ing in pseudo-trapping regions cannot complete a longitudinal 
drift and are lost from the magnetosphere (Roederer, 1967 and 
1968). 


Figure 2— Drift paths (equatorial-crossing altitudes) of electrons 
for particles mirroring near equator (ERA- 90°) and those 
mirroring at high latitudes (EPA * 0). Computations were per- 
formed assuming a magnetosphere termination at 10 R e on the 
noon meridian and a neutral sheet extending from 8 to 200 
on the midnight meridian. Tail field of 157 was assumed. 




Figure 3— Schematic representation of angular distributions on the magnetic equator for: 
(a) local midnight and (b) local noon (Roederer, 1967 and 1968). 


and the characteristics of the electron fluxes associated with varying degrees of magnetic activity. Inspection of all data available 
to us (December 6, 1966 to January 1, 1969) indicates that such sequences of quiet and disturbed electron fluxes are typical and 
can serve as a model for the discussion below. 
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Figure 4 illustrates a decrease in the energetic electron flux associated with a moderate magnetic disturbance. Typical diurnal 
variations in the electron fluxes are found on December 31 (K = 6-) and January 2 (K p = 14”). The small disturbance occurring 
on January 1 (K - 24”) created a very different curve of flux versus local time, one rich in short time-scale fluctuations and with 
energetic (E > 1.9 MeV) electron flux depressed by about an order of magnitude. 

The data for January 6, ( K p = 1 lo), presents a striking example of transition from very smooth count rates to a profile of 
count rate versus time containing a wealth of fine structure (Figure 5). Here there clearly is a pronounced diurnal variation, rang- 
ing from about a factor of two for electrons >300 keV (top trace) to a factor of 10 for electrons >1.9 MeV (bottom trace). 

The small-scale fluctuations are real and represent drift periodic echoes. Small changes in the position of the magnetopause 
result in motion of the electron flux distribution past the satellite, repeating with a period characteristic of the longitudinal drift. 

A complete analysis of this effect has been made by Brewer et al. (1969) using the data of Lanzerotti et al. (1967). An excellent 
example of drift periodic echoes for the day 351, 1966 is represented in Figure 6. 

A very disturbed period starts on January 7 (K p = 27+). Subsequent days (Figure 7) show the postdisturbance period and the 
gradual return to quiet magnetic conditions. Recoveries (and depressions) of particle fluxes appear to contain adiabatic as well as 

nonadiabatic phases. 

Examination of the data for several years obtained by this experiment shows that a more or less smooth diurnal variation per- 
sists during magnetically quiet times with the maxima in the electron flux occurring near local noon and the minima near local 
midnight. Magnetic storms introduce large changes in the electron fluxes at all local times and disrupt the diurnal variation found 
during quiet and moderately disturbed times. As viewed by the satellite, the fluxes of energetic electrons (>300 keV) are de- 
pressed and fluctuate irregularly during a storm. Another signature of a magnetic storm is a depressed and irregularly fluctuating 
magnetic field. The fluxes recover to prestorm levels with time constants that are of the order of days and appear to depend, 
among other things, on the energy of the particle. During a magnetic storm the omnidirectional electron fluxes show a great deal 
of fine temporal structure on time scales of minutes. The storm-time behavior of energetic electrons is illustrated in Figure 8. 

The unusual feature of the storm of January 13 and 14 is 
that the energetic electrons vanished for a time period of the 
order of half an hour, that is, the counting rates dropped to 
what might be expected from cosmic ray background (Figure 
9) and reappeared (Figure 10) in a manner that suggests that 
the satellite moved in and out of the region of trapping of 
energetic electrons (see compilation, 1968). This observation 
was confirmed by the UCLA magnetometer data from ATS 1 . 
The field reversed direction and became somewhat disordered 
in association with the disappearance of energetic electrons. 
This phenomenon occurred at local times from 1400 to about 
1500; previous measurements by other satellites have shown 
that the domain of energetic electrons extends to > 8 R e on 
the sunward side of the magnetosphere. 

This magnetic storm does not differ from the previous 
disturbed period (January 7-8, 1967) or from other magnetic 
storms that were observed during the first part of 1967 as far 
as surface magnetograms and ATS 1 energetic electron data are 
concerned. Typical behavior of energetic electrons during dis- 
turbed times is characterized by the following features: (a) 
transition from smooth flux-versus-time profiles to those con- 
taining much temporal structure in all integral electron chan- 
nels, (b) general depressions of energetic electron fluxes during 
the storm, with all local times affected, and (c) a poststorm re- 
covery period without any diurnal variation during which the 
electron fluxes gradually recover to the prestorm levels. The 
complete disappearance of energetic electron count rates near 
0007 UT on January 14 (Figure 9) and the subsequent appear- 
ance and disappearance of these particles (Figure 1 0) suggest 
that the satellite was moving in and out of the region of trapping. 



Figure 4— Counting rates of energetic electrons above two 
thresholds from December 31, 1966 through January 2, 1967. 
To convert from UT to local time for all ATS 1 data, subtract 

10 hours. 





Figure 5—Electron counting rates for January 6, 1967 versus 
universal time. Traces, from top to bottom, represent elec- 
trons >300 keV, >450 keV, >1.05 MeV, and >1.9 MeV. To 
convert from UT to local time for all ATS 1 data, subtract 10 

hours. 
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Figure 6— Drift periodic electron echoes observed by the Bell Telephone Laboratories 
experiment aboard ATS 1. Each repeated letter represents the same group of electrons 
making repeated transits past the satellite as they drift in longitude (Brewer et a!., 
1969). 



Figure 7— Counting rates of energetic electrons from January 
9, 1967 through January 13, 1967. This sequence of data 
illustrates the recovery of the electron fluxes after the mag- 
netic storm of January 7 and 8, 1967. 


Data Reduction 



Figure 8— Omnidirectional electron fluxes above two integral 
thresholds during the magnetic storm of January 13 and 14, 
1967. Possible crossings of the boundary of the magnetosphere 
occurred between 0008 and 0105 on January 14, 1967. 


It is clear from a consideration of the few samples of data 
discussed above, that enormous variations in electron fluxes 
are typical of the synchronous orbit. In order to prepare a representative picture of the energetic electron environment at the 
synchronous altitude, it was necessary to reduce a great body of data and present the results in a simple parametric format. 
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Figure 9— Details of the rapid electron-flux changes above two 
thresholds together with field components from the UCLA 
magnetometer data for 0000 UT to 0010 UT on January 14 , 
1967, The H{z ) component of the magnetometer data is 

shown. 



Figure 10— Details of the rapid electron-flux changes above two 
thresholds together with field components from the UCLA 
magnetometer data for 0040 UT to 0115 UT on January 14, 
1967. The H(z) component of the magnetometer data is 
shown. 


The data were reduced by computer and presented as 
plots of flux versus universal time for a given day. The aver- 
ages of the flux for each hour centered on the (UT) hour were 
also computed and ordered versus day number. An example of 
such a plot of hourly flux averages versus day number is given 
in Figure 1 1 . This figure very clearly displays the large tempo- 
ral flux variations and suggests that a model environment is 
best cast as a probabilistic model. These hourly averages are 
the basis of the average energetic-electron environment pre- 
sented in this report. Averaging tends to mask the short-term 


fluctuations in electron fluxes that are often observed, especially during magnetically disturbed periods. But since we are con- 
cerned with a long-term semiprobabilistic interpretation of this data, the treatment of data on an hourly average basis does consid- 
erably simplify data reduction. 


UT = 1000 hr 
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Figure 11 —Plot of hourly averages of omnidirectional electron (upper trace >1.05 IVIeV,and lower trace >1.9 MeV) 
for time span considered in this report. Local time of all data is midnight. Daily sum of K is plotted at the bottom 

of the figure.. 




Data were then ordered in local time and in terms of magnetic activity, as indicated by the K index. Daily sums of K p > 24 
were considered indicative of magnetic disturbance. The threshold of 2 K p > 24 was chosen because inspection of a large subset of 
the data revealed that at about this level of magnetic disturbance the diurnal variation in the electron flux levels was not visible 
above the large temporal variations. The data were analyzed in a manner first employed by Vette and Lucero (1967). We calcu- 
lated from the experimental data the probability P(F>F X ) of observing a flux greater than F x for an integral energy threshold F. 
(Figure 12). A more complete set of P(F> F ) compilations can be found in Paulikas et al. (1969). The mean of the logarithm of 
the flux 00 and the standard deviation ( o ) are calculated and the P(F> F x ) data are fitted, assuming a gaussian distribution with 
parameters jjl and a. The relevant formulas are 

number of data points with F > F 

P(F> F ) = 

x total points in sample 

i i 

° 2 = ~f ^ (log 10 F i ~ > 

i 


and 
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l exp 

r x v 

2a 2 


dZ, 



2 a 1 


The fits to the data are quite satisfactory (Figure 13), so that it is possible to characterize the electron flux distribution by the 
parameters fj. and u, which are, of course, functions of local time, electron energy, magnetic activity, and so forth. 


Local Time Variation of Electron Fluxes 

The data shown in Figures 12 and 13 summarize all the data available. We have separately computed j u and a for each local 
time. Figures 14 through 21 present plots of fi and a as functions of local time. Figure 22 summarizes the data and explicitly dis- 
plays the larger diurnal variation for the more energetic electrons (Paulikas et al., 1968). The figures in this section dearly show 
that the maxima and minima in the electron flux do not occur precisely at local noon and local midnight, respectively. At the 


DAY 340, 1966 TO DAY 60, 1968 



Figure 12— Plot of P{F> F x ) (the probability of observing an 
omnidirectional flux F greater than F x ) versus F x for four 
integral energy thresholds. 


DAY 340, 1966 TO DAY 60, 1968 



Figure 13— Fits to data presented in Figure 12. A fit has been 
made to the integral probability P{F>F x ); in addition, the 
probability density P(F) has also been plotted. 


57 




LOCAL TIME (hr) 


LOCAL TIME (hr) 


14 16 18 20 22 0 2 4 6 8 10 12 14 14 16 18 20 22 0 2 4 6 8 10 12 14 



>300 keV, versus time. Data covers Day 340, 1966 through 

Day 60, 1968. 

Figure 15— Standard deviation (a) of logarithm of electron flux 
for electrons, >300 keV, versus time. Data covers Day 340, 
1966 through Day 60, 1968. 
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Figure 16— Mean of the logarithm of the flux for electrons, Figure 17— Standard deviation (a) of logarithm of electron flux 

>450 keV, versus time. Data covers Day 340, 1966 through for electrons, >450 keV, versus time. Data covers Day 340, 

Day 60, 1968. 1966 through Day 60, 1968. 


synchronous orbit the magnetosphere is no longer symmetric about the sun-earth line, but displays about a 15° twist in the west- 
ward direction. The cause of this twist may be the draping of the interplanetary field (arriving at 1 AU at an angle of 45° to 55° 
with respect to the sun-earth line) over the magnetopause (Walters, 1964) or an effect internal to the magnetosphere caused by 
partial ring currents (Frank, private communication). Note also that a reaches a maximum near local midnight reflecting the 
greater variability of the electron fluxes on higher field lines. 
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Figure 18— Mean of the logarithm of the flux for electrons, 
>1050 keV, versus time. Data covers Day 340, 1966 through 
Day 60, 1968. 
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Figure 20-Mean of the logarithm of the flux for electrons, 
>1900 keV, versus time. Data covers Day 340, 1966 through 
Day 60, 1968. 
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Figure 19— Standard deviation (a) of logarithm of electron flux 
for electrons, >1050 keV, versus time. Data covers Day 340, 
1966 through Day 60, 1968. 
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Figure 21— Standard deviation (o) of logarithm of electron flux 
for electrons, >1900 keV, versus time. Data covers Day 340, 
1966 through Day 60, 1968. 


All data obtained by ATS 1 in the time interval of December 1966 to March 1968 are summarized in Figure 23 in the usual 
presentation of P(F > F x ) versus F x . The curves in the figure are the result of the synchronous environment compilation (AES) of 
Vette and Lucero (1967). The agreement is quite satisfactory, considering the uncertainties in the absolute geometric factors of the 
different instruments. Table 2 summarizes the means and standard deviations for the AE3 model and the ATS 1 data. 

There does not appear to be a truly major difference between the AE3 compilation (characteristic of near solar-minimum con- 
dition) and the ATS 1 data obtained primarily during 1967 (during the increasing phase of solar activity). 
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Figure 22—Summary of diurnal variation of mean of logarithm 
of flux for electrons. Data are normalized to 1100 local time. 


DAY 340, 1366 TO DAY 60, I960 



Figure 23—Comparison of ATS 1 data compilation (symbols) 
with Vette- Lucero AES mode! (solid lines). AES data for elec- 
trons, >1.9 MeV, is modified from information published in 
Vette and Lucero (1967), per private communication from 

Dr. Vette. 


Table 2 

ATS 1 , AE3 Comparison. 


Minimum 

energy 

Mean (log flux) 

Sigma (log flux) 

AE3 

ATS 1 

AE3 

ATS 1 

E e > 300 keV 

6.84 

6.795 

0.49 

0.319 

E e '>450 keV 

6.47 

6.358 

0.53 

0.396 

E e > 1.05 MeV 

5.12 

5.204 

0.615 

0.549 

E e > 1.9 MeV 

3.71 

4.138 

0.71 

0.680 


Correlation of Trapped Electron Fluxes with 
Sector Structure of Interplanetary Magnetic Field 

We have also ordered the data (for a given local time) in 
terms of solar rotations, attempting to correlate changes in the 
trapped particle population with changes in the sector struc- 
ture of the interplanetary field. Such a correlation was dem- 
onstrated for lower T-values using data from the low-altitude 
near-polar satellite 1963 38C (Williams, 1966) and was also 
presented by Vette and Lucero (1967). No such correlation is 
obvious in our data when we compare the trapped electron 
fluxes with the sector structure as reported by Wilcox and 
Colburn (1969). There are several reasons for such a lack of 
correlation. 

The rapid succession of storms during the time period of 
high solar-activity affects the particle population drastically 
through generation of new particles and loss of old particles on 
a time scale (days) that masks the long-term (about 7-day) varia- 
tion which might be expected from the effects of changing inter- 
planetary field orientation. The evolving nature of the sector 
structure during the time period under consideration also helps 
to obscure any recurrence tendencies in the electron fluxes. 

While the field geometry of low-lying {L ~ 4) drift shells 
is relatively stable, we have seen earlier that storms may totally 
alter the magnetic field from its dipolelike configuration, espe- 
cially on the nightside quadrant. Thus, the field geometry is 
not constant. 


Effect of Substorms 

The compilation of electron data presented in Figure 23 represents a group of particle fluxes which change in a correlated 
fashion; that is, an increase in the >300-keV flux is accompanied by an increase in the > 1.9-MeV electron flux. For electron en- 
ergies below about 150 keV such correlation need no longer obtain and flux changes of low-energy electrons may be anticorrelated 
with flux changes occurring among relativistic electron groups. 

An important effect at the synchronous orbit is the injection of 50 to 150 keV electrons into this region by magnetic sub- 
storms. The Minnesota group has presented a great body of results dealing with the effects of substorms on the electron popula- 
tion at the synchronous orbit (Parks et al., 1968; Amoldy and Chan, 1964; Pfitzer and Winckler, 1969; and Parks and Winckler, 
1968). The magnetospheric substorm is most often interpreted in terms of an instability which developed in the sector of dusk to 
midnight in the magnetosphere with a time scale of several hours. A general extension of the field in the anti-solar direction is 
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observed, and concurrently a loss of trapping, that is, the particle population decreases during a substorm. Recover}/ consists of a 
rapid relaxation of the magnetic field to its quiet-day configuration, with a concurrent recovery in the trapped particle fluxes. 

At this time a rapid increase of the 50- to 150-keV electron flux occurs and can be interpreted as injection of particles from 
the magnetotail. This pulse of new electrons is observed to propagate in longitude (Pfitzer and Winckler, 1969) in a manner analo- 
gous to the drift periodic echoes observed for relativistic electrons. 

During such injection of 50- to 150-keV electrons associated with a substorm, the relativistic electron fluxes do not increase 
and may even decrease. Considering the spectrum in Figure 24, one might recognize this effect as a rocking of the spectrum about 
a fulcrum at about 200 keV. The behavior of low-energy electrons thus differs in a fundamental way from that of the relativistic 
electrons. It has been suggested that the substorm injections of electrons are sufficiently frequent and powerful to provide the 
source from which the relativistic electron fluxes may be generated (Arnoldy and Chan, 1969). 

A complete description of the electron environment will thus also involve a correlation between low-energy and relativistic 
electron fluxes, that is, a probabilistic analysis of the electron spectrum. 

PENETRATION OF SOLAR PROTONS TO SYNCHRONOUS ALTITUDE 

One of the surprises contained in the initial ATS 1 observations was the fact that low-energy solar protons penetrated to the 
synchronous altitude. Cursory analysis of data obtained during the solar-proton event which began on January 28, 1967 showed 

that protons with energies greater than about 20 MeV had es- 
sentially free access to synchronous altitude at all local times. 
More detailed analysis has confirmed this result and has also 
shown that during disturbed times solar protons down to less 
than 1 MeV in energy can reach synchronous altitude without 
attenuation (Lanzerotti, 1968 and 1969 and Paulikas and 
Blake, 1969). 

Measurements of the transmission efficiency are most 
directly and unambiguously performed by comparing the flux 
outside the magnetosphere with measurements of the synchro- 
nous orbit. Because a synchronous satellite by definition cov- 
ers all local times within one 24-hour period, the transmission 
efficiency as a function of local time is easily computed. Solar 
longitude gradients in the solar-particle fluxes are assumed not 
to be important if comparisons are made with satellites which 
are outside the earth’s magnetosphere but still in earth orbit. 

A series of results which demonstrate the access of solar pro- 
tons to the synchronous orbit is presented in Figures 25 
through 28. 

The phenomena of access of solar protons to the synchro- 
nous altitude increases in complexity as one examines succes- 
sively lower energy particles. The problem is being actively 
pursued by several groups at this time, but it seems useful to 
summarize the current (fall 1969) state of knowledge for the 
purpose of this report. 

Protons with energies ^20 MeV reach, synchronous alti- 
tude without significant attenuation when compared with 
fluxes above the same threshold outside the magnetosphere. 
Thus, observation of protons at this energy, carried out at the 
synchronous altitude are equivalent to measurements of the 
interplanetary flux and it is possible to derive such results as 
source spectra, source strengths, and interplanetary diffusion 
coefficients from these data. This result adds energetic solar- 
proton monitoring to the already numerous useful tasks which 
can be performed by satellites at the synchronous orbit. 
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ELECTRON ENERGY (MeV) 


Figure 24-Spectrum of electrons at synchronous altitude ob- 
tained aboard OV2 5 (J. R. Stevens and A. L. Vampola, private 
communication). Substorms primarily affect particle fluxes 
with energies below 200 keV. 






UT (hr) 


Figure 25— Comparison of ATS 1 and IMP 3 data obtained 
during initial phase of January 28, 1967 solar proton event. 
IMP 3 directional fluxes have been converted to omnidirec- 
tional fluxes assuming isotropy. The assumption of isotropy is 
reasonable since the source region of particles was on back side 
of the sun. 



JANUARY 29 JANUARY 30 JANUARY 31 FEBRUARY I 

(DAY 29) (DAY 30) (DAY 31) (DAY 32) 


Figure 26— Transmission coefficient for 5- to 21-MeV and 21- to 70-MeV protons for 
January 29 through February 1, 1967. Uncertainties in geometric factors are the main 
source of error in computations of transmission efficiencies and amount to ±30% for 
the 5- to 21-MeV channel and ±50% for the 21- to 70-MeV channel. Also shown 
are the 3-hour K indices and hourly DST values (M. Sugiura and S. Cain, private 
communication). 







Figure 27— Comparison of low-energy solar protons {dots) and alphas (triangles) measured inside the 
magnetosphere aboard ATS 1 (open symbols) and aboard Explorer 34 outside the magnetosphere 
(solid symbols) by the BTL experimenters (Lanzerotti, 1968). 


Protons with energies below about 20 MeV are affected by the magnetosphere. This interaction manifests itself as a diurnal 
variation in the solar-proton flux caused by the asymmetric nature of the magnetosphere. The solar-proton diurnal variation is in- 
verted with respect to the trapped electron variation. More solar protons are observed near local midnight than near local noon. 
During magnetic quiet and near local noon, the flux at the synchronous altitude of 5- to 21 -MeV protons is decreased by a factor 
of five to ten with respect to interplanetary levels. However, at night or during disturbed times the attenuation is less than a factor 
of two. Direct and indirect evidence indicates that solar protons below about 20 MeV are temporarily trapped in the geomagnetic 
field. The trapping lifetime is of the order of half an hour or less. Because the particles are trapped, they may be locally acceler- 
ated by changes in the geomagnetic field, and the flux (in a given energy interval) may temporarily increase to exceed the inter- 
planetary flux (Figure 29). 

Protons of a few MeV energy are still able to reach the synchronous altitude. The ease of access appears to depend in a com- 
plicated way upon particle energy and the state of agitation of the magnetosphere. Occasionally, significant temporal delays 
(hours) are observed between arrival of solar particles in interplanetary space near the earth and arrival of those particles at the 
synchronous altitude (Lanzerotti, 1969a). 

The behavior of various energy groups of solar protons during the event of January 28, 1 967 is illustrated in Figure 30. Note 
that diurnal variation is most pronounced for particles of intermediate energy; the diurnal variation is absent for protons E > 20 
MeV and also for protons with E p =€ 3 MeV. 
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Figure 28— Transmission coefficient for 5- to 21-MeV protons 
as a function of 3-hour K . Data in the 0900 to 1500 local 
time quadrant (day) and the 2100 to 0300 local time quadrant 

(night) are separately identified. 



2700 28.00 29.00 30.00 31.00 32.00 33.00 


ur (hr) 

Figure 30— Behavior of solar proton fluxes in five energy inter- 
vals at the synchronous altitude. Data were obtained by the 
BTL experiment aboard ATS 1 during the solar proton event 
of January 28, 1967. Factor in parentheses is scale factor 
appropriate for each plot (Roberts et al., 1968). 
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Figure 29— Effect of sudden commencement on solar proton 
fluxes. Solar proton flux on July 13, 1968 (Day 195). Shown 
are the ATS 1 proton flux 10-minute averages in two energy 
intervals (5- to 21-MeV and 21- to 70-MeV) and the Explorer 
34 (> 1 0 MeV) counting rate. (Explorer 34 was outside the 
magnetosphere.) Aiso shown (as an indicator of the state of 
the magnetic field fluctuations) is the trapped flux of electrons 
>1.05 MeV. The proton flux enhancement following the sud- 
den commencement is followed by a flux increase associated 
with energetic storm particles (Paulikas and Blake, 1969). 

It is important to note that orbit tracing in realistic, but 
static, models of the magnetosphere does not predict the ob- 
served situation regarding solar-proton access to the synchro- 
nous altitude (Lanzerotti, 1969a and Paulikas and Blake, 1969). 
Only by making extreme assumptions regarding the distance to 
the magnetopause or the strength of the tail field in, for exam- 
ple, a Williams-Mead model can protons of about 1 0 MeV mir- 
roring near the equator reach 6.6 R e according to such calcula- 
tions. It is necessary to invoke a diffusion mechanism, as yet 
unspecified, which transports particles radially deeper into the 
magnetosphere than predicted by orbit tracing alone. Direct 
diffusion into the flanks of the magnetosphere has been pro- 
posed by Lanzerotti (1968). An alternative mode of rapid 
access can be visualized by noting that (Figure 1) the quasi- 
trapped regions on the night side, at somewhat larger geocentric 
distances than the synchronous orbit, are accessible from the 
flanks of the magnetosphere by drift alone. If pitch-angle scat- 
tering of the particles is added, radial diffusion will occur be- 
cause of the dependence of drift path on pitch angle. 

Solar-proton observations aboard ATS 1 during 1967 are 
summarized in Figure 3 1 . The data have been ordered in terms 
of solar rotations. Note that solar proton fluxes above back- 
ground are present a significant fraction of the time; however, 
the major contributions to the integrated flux for 1967 were 
made by the events of January 28 and May 23. 
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Figure 31-Compilation of (5- to 21-MeV) solar proton fluxes observed aboard ATS 1 during 1967 . 
Data are ordered in terms of solar rotations (Paulikas and Blake, 1968). 









SUMMARY 


This paper has presented a state-of-knowledge report on the trapped energetic electrons and energetic solar-protons observed 
at the synchronous altitude. Important progress has been made in understanding the behavior of the particle population at this 
altitude, although large gaps still exist in our knowledge. The increase in knowledge has been exceptionally rapid because of the 
availability of a large volume of high-quality (nearly continuous) data from ATS 1 and other satellites. An important ingredient in 
this progress is the relative ease of data analysis for experiments operating at the synchronous altitude. 

There exists still only fragmentary data regarding pitch-angle distributions of energetic particles. Interesting new results can 
be expected from analysis of pitch-angle data. 

Some progress is also being made in quantitative understanding of the particle-field interactions which operate in this region 
of the magnetosphere. It is this kind of work that promises to unravel the fundamental processes in the radiation belt and build a 
sound theoretical basis for the phenomena presented in this report. 
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